The 

Computation of 
Fixed Points and Applications 



Michael J. Todd 



124 



Springer 




Lecture Notes 
in Economics and 
Mathematical Systems 

Managing Editors: M. Beckmann and H. P. Kunzi 
Mathematical Economics 



124 



Michael J. Todd 



The Computation of Fixed Points 
and Applications 




Springer-Verlag 
Berlin Heidelberg GmbH 




Editorial Board 

H. Albach • A. V. Balakrishnan • M. Beckmann (Managing Editor) 
P. Dhrymes • J. Green • W. Hildenbrand • W. Krelle 
H. P. Kunzi (Managing Editor) • K. Ritter • R. Sato • H. Schelbert 
P. Schonfeld 

Managing Editors 
Prof. Dr. M. Beckmann 
Brown University 
Providence, Rl 02912/USA 

Author 

Michael J. Todd 

School of Operations Research 
and Industrial Engineering 
Cornell University 
Ithaca, New York 14853/USA 



Library of Congress Cataloging in Publication Data 

Todd, Michael J 1947- 

The computation of fixed points and applications. 

(Lecture notes in economics and mathematical sys- 
tems ; 124) (Mathematical economics) 

Bibliography: p. 

Includes index. 

1. Fixed point theorems (Topology) 2. Triangu- 
lating manifolds. 3,. Economics, Mathematical. 

I. Title. H. Series. III. Series: Mathematical 

economics . 

QA611.7.T63 514'. 3 76-9786 

AMS Subject Classifications (1970): 05-01,55-04, 55C20, 57C15.90A15, 
90C30 



ISBN 978-3-540-07685-8 ISBN 978-3-642-50327-6 (eBook) 

DOI 10.1007/978-3-642-50327-6 

This w.ork is subject to copyright. All rights are reserved, whether the whole 
or part of the material is concerned, specifically those of translation, re- 
printing, re-use of illustrations, broadcasting, reproduction by photocopying 
machine or similar means, and storage in data banks. 

Under § 54 of the German Copyright Law where copies are made for other 
than private use, a fee is payable to the publisher, the amount of the fee to be 
determined by agreement with the publisher. 

© by Springer-Verlag Berlin Heidelberg 1976 

Urspriinglich erschienen bei Springer-Verlag Berlin Heidelberg New York 1976. 



Prof. Dr. H. P. Kunzi 
Universitat Zurich 
8090 Zurich/Schweiz 





TABLE OF CONTENTS 



Preface v 

Notation v j_ 

Chapter I Brouwer’s Theorem 1 

Chapter II Some Applications of Brouwer’s Theorem 17 

Chapter III Triangulations 24 

Chapter IV Algorithms to Find Completely-Labelled Simplices 40 

Chapter V Extensions of Brouwer's Theorem 54 

Chapter VI Applications of Kakutani’s Theorem and its Extensions .... 64 

Chapter VII Eaves' First Algorithm 73 

Chapter VIII Merrill’s Algorithm 82 

Chapter IX Homotopy Algorithms 94 

Chapter X Triangulations with Continuous Refinement of Grid Size .... 101 

Chapter XI Measures of Efficiency for Triangulations 116 



References 



126 




PREFACE 



Fixed-point algorithms have diverse applications in economics , optimization , 
game theory and the numerical solution of boundary -value problems. Since Scarf’s 
pioneering work [56,57] on obtaining approximate fixed points of continuous mappings, 
a great deal of research has been done in extending the applicability and improving 
the efficiency of fixed-point methods. Much of this work is available only in 
research papers, although Scarf’s book [58] gives a remarkably clear exposition of 
the power of fixed-point methods. However, the algorithms described by Scarf have 
been superseded by the more sophisticated restart and homotopy techniques of Merrill 
[48,49] and Eaves and Saigal [14,16]. To understand the more efficient algorithms 
one must become familiar with the notions of triangulation and simplicial approxima- 
tion, whereas Scarf stresses the concept of primitive set. 

These notes are intended to introduce to a wider audience the most recent 
fixed-point methods and their applications . Our approach is therefore via triangula- 
tions . For this reason , Scarf is cited less in this manuscript than his contribu- 
tions would otherwise warrant. We have also confined our treatment of applications 
to the computation of economic equilibria and the solution of optimization problems. 
Hansen and Koopmans [28] apply fixed-point methods to the computation of an invariant 
optimal capital stock in an economic growth model. Applications to game theory are 
discussed in Scarf [56,58], Shapley [59], and Garcia, Lemke and Luethi [24]. 

Allgower [1] and Jeppson [31] use fixed-point algorithms to find many solutions to 
boundary-value problems. Infinite-dimensional cases are also discussed by 
Freidenfelds [21] and Wilmuth [73]. The Schauder Projection theorem (see 
Freidenfelds [21] or Smart [61]) describes how a finite-dimensional approximation 
can be obtained. Our treatment is confined to finite-dimensional spaces throughout. 
More recent developments we have been unable to cover are the orientation theories 
of Shapley [60], Lemke and Grotzinger [44], Todd [66] and most generally Eaves [15] 
and Eaves and Scarf [17]; and the algorithmic improvements and convergence analysis 
of Saigal [53,54]. Eaves [15] contains a comprehensive bibliography. 

The manuscript is organized as follows. Chapter I gives a classical (non- 
algorithmic) proof of Brouwer's theorem from Sperner's lemma, thus introducing the 
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reader to several important concepts. A number of applications are described in 
Chapter II. We provide a formal treatment of triangulations in Chapter III with 
descriptions of some important particular triangulations. The latter are used in 
Chapter IV in algorithms for computing approximate fixed points of continuous func- 
tions. The applications of Chapter II motivate extensions from functions to point- 
to-set mappings. Chapters V and VI parallel Chapters I and II in proving and 
applying Kakutani's fixed-point theorem. Chapter VII describes an algorithm for 
computing Kakutani fixed points. This algorithm and those of Chapter IV are ineffi- 
cient if a good approximation is desired. In Chapters VIII and IX we describe the 
more sophisticated restart and homotopy algorithms. The latter require special 
triangulations which are developed in Chapter X. Finally, Chapter XI describes some 
measures that can be used to compare different triangulations when used for fixed- 
point computation. 

We have included a number of challenging exercises to increase the reader’s 
understanding of the material. Some numerical examples of the algorithms are given 
in the text. The reader is assumed to be familiar with real analysis and linear 
programming, including lexicographic resolution of degeneracy. We also assume the 
Kuhn-Tucker conditions for nonlinear programming known. 

This manuscript arose out of a course in computing fixed points that I gave at 
Cornell University in spring 1975. I am grateful to Michel Cosnard, Pierre Dejax, 
Pradeep Dubey, Etienne Loute, Shigeo Muto, Bob Rovinsky and Prakash Shenoy for 
preparing excellent notes. The National Science Foundation, through grant GK-42092, 
provided support during the preparation of this manuscript. I would like to thank 
Mrs. Kathy King for her excellent typing. Finally, my thanks go to my wife Marina 
for her encouragement and assistance. 

Notation 

N: Set of integers {l,2,...,n}. 

Nq'. Set of integers {0,l,...,n}. 

R m : Set of m-dimensional real column vectors, with coordinates generally indexed 

1 through m. However, the coordinates of R n+ ^" are indexed by N^. 
u 1 : ith unit vector in R n , i £ N; u = u 1 . 
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jth unit vector in R 1 



n+1 



j € N 



O’ 



y. v 3. 



R^: Nonnegative orthant of R m ; {x € R m |x >_ 0}. 

CUD, COD, C^D: Union, intersection and set difference of the sets C and 

C + D, C - D: Algebraic sum and difference of sets C and D in R™ ; 

{c + d|c £ C, d £ D} and {c - d |c £ C, d £ D} respectively. 



D. 



AC: {Ac|c £ C}. 

Ilxjl^: Euclidean norm of the vector x £ R m ; x?)^"^. 

li x lL : ^oo norm x £ R m ; max.Jx.jJ. 

||A ||p : The operator norm; max{||Ax||p | ||x ||^ = 1} for p = 2,® 

real kxm matrix. 



B m : Euclidean unit ball in R m ; {x £ R m | ||x < 1}. 




B(x,p) : 


Ball with center x radius p; {x} + pB m if x 


€R m . 


B(C,p ) : 


C + pB m if C C R m . 




C: Closure of C; f) {B(C,e)|e > 0}. 




int C: 


Interior of C; {x £ c|3 e > 0 with B(x,e) C C}. 




diam C: 
P 


Diameter of C; sup{ ||x-y ||Jx,y £ C}, p = 2 or 


°°* 


mesh G: 
P 


Mesh of G; supidiam^ClC £ G} for p = 2 or ®. 


, where 



where A is a 



G is a family 



of subsets of R m . 




CHAPTER I: BROUWER’S THEOREM 



I .1 . Probably the most famous fixed-point theorem states that a continuous function 
from a closed n-cell to itself leaves at least one point fixed. The Dutch mathema- 
tician L. E. J. Brouwer proved this result in 1912 [5] using degree theory. An equiv- 
alent theorem in the case of differentiable functions was proved earlier by Bohl [4], 
who used Green's theorem. Our concern is more with the computation of fixed points 
than their existence; we will follow a later approach based on the purely combina- 
torial lemma of Sperner [62]. This approach is closest to the algorithms we shall 
develop, and the machinery will be of value later. However, to avoid some of the 
cumbersome details, we will give only an intuitive idea of the notions of simplex and 
triangulation. A rigorous treatment appears in Chapter III. 

In this section we state Brouwer's theorem formally. Section 2 shows that a 
proof for the standard simplex is sufficient and gives some examples suggesting var- 
ious methods of proof. In Section 3 we reduce Brouwer's theorem to Sperner 's lemma, 
which is proved in Section 4. 

1.1 Definition. A function h is a homeomorphism if it is one-to-one and onto, 
and both h and h ^ are continuous. A closed n-cell is a homeomorphic image of 
B n , i.e., C is a closed n-cell if there is a homeomorphism h: B n C. 

1.2 Theorem (Brouwer). Let C be a closed n-cell and let f: C ** C be contin- 
uous. Then f has a fixed point, i.e., there is an x* £ C with f(x*) = x*. 



1.2 . In this section we show that is is sufficient to prove 1.2 when C is the 
standard simplex and give some examples of the necessity of the conditions and possi- 
ble methods of proof. 



2.1 Definition . The standard simplex S n is the convex hull of v^v 1 ,...,^ 

in R n , i.e., S n = {x £ R^ +1 |v T x = 1}. For i £ N Q , S? denotes the face of S n 

opposite v 1 , i.e., {x £ S n |x. = 0}, and the boundary of S n is 9S n = U. M S?. 

1 0 
We show below that S D is a closed n-cell, but to get an intuitive feel for 



this concept , we first give an easily visualized subclass : 
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2 . 2 Lemma . If C C R n is compact and convex with a nonempty interior, then 
C is a closed n-cell. 



Proof . We construct a function h: B n -*■ C and show that it is one-to-one and 
onto. The proof that h and h" 1 are continuous is omitted. 

Pick c £ int C. For 0 ^ d ^ R n define 0(d) = max{0 £ R|c + 0d £ C}. (The 

maximum is attained because C is compact.) Since c £ int C, 0(d) > 0, and 
0(Ad) = A ^0(d) for A > 0. (0 is closely related to the Minkowski functional; 
see [61].) 

Now define h: B n -► C by setting h(d) = c + ||d|| 2 0(d)d if d ^ 0, and 

h(0) - c. Clearly, h(d) = c iff d = 0, so that h(d) = h(d’) = c => d = d' . Now 

let h(d) = h(d’) i c; then ||d || 2 ©(d)d = ||d* || 2 0(d' )d f so that d’ = Ad, A > 0. 

But then ||d|| 2 0(d)d = A||d|| 2 A _1 0(d)Ad, and A = 1. Thus h is one-to-one. 

We now show that h is onto. Since h(0) = c, we only need to find d £ B n 
so that h(d) = x for x £ C, x i c. Let d = (x-c)/||x-c|| 2 0(x-c) . It is easily 

checked that d £ B n and h(d) = x. 

We use 2.2 to show that S n is a closed n-cell by introducing the set C n = 

(x £ R n 1 1 >_ x^ >_ . . . >_ x^ >_ 0}, which will appear in Chapter III. 

2.3 Lemma. S n is a closed n-cell. 



Proof . Since by 2.2 C is a closed n-cell, we need only show that S and 
C n are homeomorphic . Define h: C n -► S n and h S n -> C n by h(c) = + Qc, 



h (s) = Q t s with Q the (n+l)xn matrix 

~ 0 1 ... 1 ^ 

matrix 



• o 

+i • 

0 • 



-1 



and Q f the nx(n+l) 



Obviously, h is one-to-one and onto, and h and h 

* • • 

0 ... 0 1 

are both continuous. Thus S n is homeomorphic to C n and hence to B n . Q] 



-1 



The following result shows that it is sufficient to prove Brouwer’s theorem for 
S n . While this simplifies the proof, it is certainly not helpful when one wants to 
compute a fixed point of some other closed n-cell, as one needs explicitly a 
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homeomorphism between the two. This problem will be satisfactorily dealt with when 
we treat Kakutani's theorem in Chapter V. 

2 . 4 Lemma . Brouwer’s theorem is true for all closed n-cells C if it is true 
for S n . 



Proof . Let f: C C be continuous. We must show that f has a fixed point. 
Since C and S n are closed n-cells, we have homeomorphisms h: B n C and 
h Q : B n -*■ S n . Then the composite function f Q = h Q h ^fhh^ 1 : S n -► S n is continuous 
and has a fixed point x* by hypothesis. Hence fChh^Cx*)) = hh^Cx*) and f has 
a fixed point. | | 

2 . 5 Examples of the necessity of the conditions of 1.2. 

(i) C is not a closed n-cell (see also exercise 5.1). 

2 2 

(a) C is closed and int C i 0, but C is not convex. Let C = 2B B 
and f(x) = -x. Then f rotates the torus by tt and is clearly continuous, but 
has no fixed points. 

f ( x ) 



2 

(b) C is convex and int C t 0, but C is not closed. Let C = int B and 
f(x) = ^ x + y u 1 . Clearly f is continuous and has no fixed points. 
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(ii) f is not continuous. 
Let C = B 2 and f(x) 



0 if x 1 0 

. Then f has no fixed points. 

u 1 if x = 0 



2.6 Possible Methods of Proof . 

(a) (n = 1) Since S 1 is homeomorphic to [0,1], consider a continuous func- 
tion f : [0,1] -► [0,1]. If f(0) = 0 or f(l) = 1, we are done. Otherwise, let 
g(x) = f(x) - x; then g is continuous with g(0) > 0 > g(l). By the intermediate 
value theorem, g has a zero x* in [0,1] and hence x* is a fixed point of f. 
Intuitively, the graph of f must cross the diagonal of the unit square, giving a 




It is not clear how the argument can be extended to higher dimensions — but see 

(c). 

2 2 

(b) (n = 2) Let us assume we have a continuous function f: B -*■ B without 
fixed points and try to establish a contradiction. We can then define a function 
h: B 2 -► 9B 2 E {x £ R 2 |||x|| 2 = 1} by setting h(x) to be the point of 3B 2 on the 
line from f(x) to x: 
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If f has no fixed points, it can be shown that h is continuous, and clearly h 

2 

leaves fixed each point of 8B . It seems intuitively clear that no function can 
2 2 

carry B into 9B , leaving the latter fixed, without "ripping" the interior of 
2 

B and thus losing its continuity. One can make this argument rigorous for any n, 

using either homology or the more elementary proof of Hirsch [30]. Hirsch in effect 
-1 2 

traces the inverse image h (b) of a point b £ 3B and shows that it can only 
2 

"disappear" inside B . It must disappear near a fixed point of f; Kellogg, Li, 
and Yorke [35] have recently constructed an algorithm based on this idea — they require 
f to be twice continuously differentiable. We will not pursue this approach. 

(c) We now motivate the method we will follow. Consider the case n = 1. Let 
f: S 1 be continuous. For i = 0,1, let be the set of points of 

whose i^ coordinate does not increase under f. Clearly, each point of lies in 

either C Q or C^; v° £ C^ and v 1 £ C^; and C Q and C^ are closed. If C Q 
and C^ intersect, then any point in their intersection must be a fixed point of f, 
for neither of its coordinates can increase under f while their sum remains 1. 

Thus the first step in our proof of Brouwer’s theorem is to reduce it to the lemma 
of Knaster, Kuratowski and Mazurkiewicz , stating that under certain conditions a 
family of sets must have a nonempty intersection. Unfortunately, this lemma seems 
no easier to prove than the original theorem, and we will have to reduce it further 
to a purely combinatorial lemma. Consider again the case n = 1. If we can find 
arbitrarily close pairs of points, one in C^ and one in C^, then since C^ and 
C^ are closed and is compact, our lemma will be established. To provide an 

abundant supply of pairs of points, we can divide into a family of small seg- 

ments . An endpoint of a segment will be labelled 0 or 1 according to whether it lies 
in Cq or C^. (If it lies in both, the lemma is proved.) Since one end of S^ 
has the label 0 and the other the label 1, it is clear that there is a small segment 
whose endpoints are labelled 0,1, thus yielding a close pair of points. We have 
the following picture : 
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In higher dimensions we must find (n+l)-tuples of close points, one in each C^. 

Again we obtain them by dividing S n into small pieces, called simplices. We divide 
2 3 

S into triangles, S into tetrahedra, and so on. In this way we reduce our lemma 
about families of closed sets to Sperner’s lemma concerning labelled simplicial sub- 
divisions . 



1.3. In this section we reduce Brouwer's theorem for S to Sperner's lemma. Al- 
though this step can be made directly (3.7), we proceed more slowly (and, we hope, 
intuitively) via the Knaster-Kuratowski-Mazurkiewicz (K-K-M) lemma [36]. (This lemma 
has other interesting applications; see exercise 5.3.) 



3 . 1 Lemma (K-K-M). Let C^, i £ N Q , be a family of closed subsets of S 
satisfying the following conditions: 

(i) S n = U.,„ C.; and 

i£N 0 1 

(ii) If MICH 0 and J = N 0 ^ I, then H. x s'? c U j6J C.. 

Then 0 C. jf 0. 
i€N 0 i 

For n = 1, (ii) says only that and ; for n = 2, the con- 

ditions are illustrated by the picture below. The shaded section shows that 



n. „ c. 0 0. 

* N 0 1 
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Proof . Let f: S n S n be continuous. For i £ N Q let 

c! = {x fS n |f.(x) < x. > 0} and C. = C.. We now show that the C., i f N_, sat- 

i 1 l — l li i ^ 0 

isfy the hypotheses of the K-K-M lemma. If x £ S U lies in none of the C^, i £ N Q , 

then f.(x) > x. for i £ I = {i £ N. lx. > 0}. Then 1 = v^x = Y T x. < Y_ f.(x) < 
l i O' l L Ii L Ii — 

v T f(x) = 1, a contradiction. Hence S n c U. „ c! c [L XT C.. Since for i c N_, 

- ^Nq i - i£N Q i c O’ 

CS n , condition (i) is verified. Also, by definition, if x £ S*?, x c|, 

i € I. Thus for J = N_ ^ I, x £ U. . c! C U. T C., and condition (ii) is verified. 

0 1 ~ 3 

If the K-K-M lemma holds, we can deduce the existence of x* f fl. .. C.. Since 

i£H 0 1 

f is continuous and x* f C ! , i f N_, we have f.(x*) < x*, i f N. Now 1 = 

i 0 i — i 0 

v T x* = v^f(x*) shows that f.(x*) = x*, i £ N_, and so x* is a fixed point of 

i l 0 r 

□ 



Note: Scarf [58] has proved a lemma similar to 3.1 in which condition (ii) is 

replaced by: 

(ii)' For all i £ N A , s"cC.. 

0 l—i 

One can easily verify that this version also implies Brouwer’s theorem for S n . 
Freidenfelds [22] has shown that the two versions are equivalent. 



3.3 Simplices and Triangulations . As we discussed in 2.6(c), we shall prove 
3.1 by exhibiting, for any £ > 0, an (n+1) -tuple of points of S n within e of 
one another, with one point in each C^. We produce these (n+1) -tuples by dividing 
S n into many small simplices. The concepts of simplex and triangulation will be 
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introduced informally now; we postpone rigorous definitions until Chapter III. 

A closed j -simplex is the convex hull of j+1 points in general position (af- 
finely independent) in R ; the points are called the vertices of the simplex. An 
open j -simplex is the relative interior of a closed j -simplex. 

n Closed n-simplex Open n-simplex 

0 

1 

j A 

Note: S n is a closed n-simplex. 

A face of a simplex a is a simplex all of whose vertices are vertices of a. 
Thus S? is a closed face of S n for each i £ N Q . 

Two simplices are incident if one is a face of the other. Two j-simplices are 
adjacent if they share a ( j-l)-simplex as a face. 

A triangulation G of S n is a finite collection of open n-simplices such 
that the open n-simplices, together with all their open faces, form a partition of 
S n , i.e., S n is their disjoint union. As we shall see in Chapter III, this rather 
non-intuitive definition, when couched in terms of closed n-simplices, is equivalent 
to the conditions: 

(i) the closed n-simplices cover S n ; and 

(ii) if two closed n-simplices meet, their intersection is a common face. 



(point) • 

(line segment) >— ( 

( triangle ) 



A 



( tetrahedron ) 



y lN 



Thus the following configurations are ruled out: 
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Exercise 5.5 shows that the construction of regular triangulations for n > 2 is not 
as easy as it appears for n = 1 and n = 2. 

We call a vertex of a simplex of a triangulation G a vertex of G. Hence- 
forth, a simplex will mean an open simplex, but little confusion will result from 
thinking of closed simplices. 



3.4 Properties of Triangulations of S n . The following intuitive results will 
be proved in Chapter III. 

(a) If G is a triangulation of S n and i is an (n-l)-simplex that is a 
face of an n-simplex of G, then either 

(i) t C 3S n and x is a face of just one a £ G; or 

(ii) t £ 3S n and t is a face of precisely two simplices in G. 

(b) There exist triangulations of S n of arbitrarily small mesh. (Recall that 
the mesh of G is sup diam a.) 

(c) Let G be a triangulation of S n and i £ N . Then the collection G' 

of (n-l)-simplices that are faces of simplices of G and lie in is a triangu- 

lation of S? (defined in the obvious way). 

We can now state 



3.5 Lemma (Sperner [62]). Let G be a triangulation of S n with each vertex 
of G labelled with an integer in such that no vertex in is labelled i. 

(Such a labelling is called admissible . ) Then there is a simplex in G whose 
vertices carry all the labels in N Q (a completely-labelled simplex). 



Example for n = 2. 
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The statement above is the weak form. The strong form of the lemma asserts the 
existence of an odd number of completely-labelled simplices. (There is even a 
’’super-strong” version. In the example above, note that there is one more small 

triangle with the labels 0, 1, 2 appearing counter-clockwise (the same direction as 
0 1 2 2 

v , v and v in S ) then in the reverse direction. The super-strong version 
asserts this in general, for any n, but the statement and its proof require orien- 
tation arguments.) 

For related combinatorial results, see Kuhn [37], Tucker [69], whose lemma 
relates to antipodal fixed-point theorems, and Ky Fan [18], who synthesized Sperner' s 
and Tucker’s lemmas in very general results. 

3.6 Sperner ’s Lemma (weak form) implies the K-K-M Lemma . 

Proof . Let CL , i £ N Q , be closed sets satisfying the hypotheses of the 
K-K-M lemma. Using 3.4(b), let G^, k = 1,2,..., be a sequence of triangulations 
of S n with mesh -► 0. For each k, label each vertex y of G with i = 
min{j £ | y £ CL , y £ S^}. (The existence of i follows from the conditions 

(i)-(ii) of 3.1.) This is an admissible labelling; if Sperner ’s lemma is true, we 
can deduce the existence of a completely-labelled simplex in G^. Let the 

vertices of be y^, i £ N Q , with y^ labelled i, so that y^ £ CL, 

1 € N q . 

Now the sequence y^, k = 1,2,..., lies in the compact set S n , and there- 
fore there is a convergent subsequence. Without loss in generality, assume 
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yk° + x* £ S n . Since mesh G^ -► 0, we have lim^^ y ^ = x* for each i £ N Q . 

Since C. is closed, we deduce that x* £ C. for all i £ N., which establishes 
1 i 0 

the K-K-M lemma. Q 



One might conclude that when Sperner's lemma is used to prove Brouwer's theorem, 
a completely-labelled simplex gives an approximate fixed point. The result below 
makes this precise. It is important to note that we use the term "approximate fixed 
point" to mean a point that is close to its image, while not necessarily close to a 
fixed point. 



f 

y 

G 



3 . 7 Lemma . Let G be a triangulation of 
S n S n be such that Hx-z^ < 6 implies 
of G i = min{j|fj(x) <_ x^ > Then if 

and x* £ a, ||f(x*) - x*^ <_ n(e+6). 



S n of mesh^ at most <5 . Let 
||f(x) - fCz)^ <_ e. Label a vertex 
a is a completely-labelled simplex of 



Proof. Let c have vertices y 1 , i £ N^, where y 1 has label i. Then for 



each i £ N Q we have 



f.(x*) 



X? = (f.(x*) - f . (y ) ) + (f . (y ) 
i i i i 



y-) 



(y- 



The hypotheses of the lemma guarantee that the first term on the right-hand side is 
at most e, while the last is at most 6. Since y 1 has label i, the middle 
term is nonpositive. Hence for each i £ N , 

f.(x*) - xV < e + 5. 

i l — 

T T 

Also, since v f(x*) = v x* = 1, we have for each i £ 

f.(x*) - X* = -y.,. (f.(x*) - X?) > -n(e+6). 
i i 3 3 ~ 

The conclusion now follows. | [ 



Exercise 5.6 asks the reader to develop an admissible labelling scheme that 
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allows the factor n above to be eliminated. If enough is known about the function 
f, 3.7 allows one to find a point within y of its image. Let e = y/2n and 
choose 6 so that 5 and e satisfy the condition in 3.7. Then triangulate S n 
with a triangulation of mesh at most min{<5,e} and find a completely-labelled sim- 
plex (by exhaustive search or by one of the algorithms of Chapter IV). 

1.4 Proof of Spemer's Lemma . See Lyusternik [41] for an intuitive proof for 
n = 2. We will prove the strong form of Sperner’s lemma by induction on n. We 
could start with the case n = 0, which is trivial. However, we start with n = 1; 
we also forgo some easy proofs of this case and instead use the same line of rea- 
soning we employ for the inductive step. 

4.1 The Case n = 1 . The picture is as follows: 

0 0 10 111 



Consider the incidences of vertices ( (n-l)-simplices) labelled 0 with segments 
(n-simplices ) of the triangulation. We will count these incidences in two ways. 

(i) Add up the incidences for each segment with at least one endpoint labelled 
0. Such segments split into two sets: 

A - the set of segments with one endpoint labelled 0 and the other 1; and 

B - the set of segments with both endpoints labelled 0. 

Each segment in A contributes 1 to the count of incidences, and each segment 
in B contributes 2 . Hence the total count is | A | + 2 J B [ . 

(ii) Add up the incidences for each vertex labelled 0. Such vertices also 
split into two sets : 

C - the set consisting of just v^, the only vertex in 3S labelled 0; and 

D - the set of internal vertices (not in 3S^) labelled 0. 

Using 3.4(a) (trivial in this case), we see that the vertex in C contributes 
1 while each vertex in D contributes 2 to the count of incidences . So the total 
count is |c|+2|d|=1+2|d|. 

From |a|+2|b|=1+2|d| we deduce that | A | is odd, which is precisely 
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the claim of Sperner's lemma for n = 1. 

4.2 The Inductive Step . Assume that the strong form of the lemma is true for 
dimension n-1, and let G be a triangulation of S n with the vertices of G 
admissibly labelled. Let H be the collection of (n-l)-simplices that are faces of 
simplices in G and whose vertices carry all the labels 0,1,..., n-1. We will again 
count the incidences of G and H. 

(i) Add up the incidences for each a £ G whose vertices have the labels 
0,1,..., n-1. Such simplices fall into two sets: 

A - the set of completely-labelled simplices of G; and 

B - the set of simplices of G whose vertices carry all the labels 

0,1,..., n-1 but not n. 

Each simplex in A has just one face in H, while each simplex in B has 
precisely one duplicated label among its vertices, and hence two faces in H — namely, 
those faces opposite the two vertices sharing the same label. Thus the total count 
of incidences is | A | + 2 | B | . 

(ii) Add up the incidences for each x £ H. H is partitioned into two sets: 

C - the set of x £ H with x c 9S n ; and 

D - the set of x £ H with x £ 9S n . 

By the rules of an admissible labelling, each x £ C must lie in S^. 

By 3.4(a), each x £ C contributes 1 and each x £ D contributes 2 to the 
count of incidences of G and H. Thus the total count is |c| + 2 | D | and we 
obtain | A | + 2 | B | = |c| + 2 | D | . 

To show that | A | is odd, we must show that |c| is odd. But consider the 
collection G' of (n-l)-simplices that are faces of simplices of G and lie in 
S^. By 3.4(c), G' triangulates S^. Moreover, is clearly homeomorphic to 

S n \ by the deletion or addition of the last zero coordinate. Considered as a 
triangulation of S n \ the vertices of G* are admissibly labelled, and the induc- 
tive hypothesis allows us to claim that there is an odd number of completely-labelled 
simplices in G T (i.e., simplices whose vertices carry all the labels 0 ,1 , . . . ,n-l). 
But these simplices are precisely those belonging to C; hence |c| and | A | are 
odd , and the inductive step is complete . [[[] 
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From 2.4, 3.2, and 3.6, we have finally proved Brouwer’s theorem. 



2 

4 . 3 Example (n = 2). We show below a labelled triangulation of S and iden- 
tify the sets A, B, C, and D and their incidences. 

2 

v 



2 




2-simplices in A 



Incident 1-simplices in C U D 



2-simplices in B 
a , 



Total incidences = I A I + 2 B = 9 



1-simplices in C 

T i 

T 2 

T 3 

1-simplices in D 

\ 

T 5 

T 6 



Incident 2-simplices in A U B 



a 

o 

a 



1 * 

2 ’ 

3’ 



a 

a 

a 



2 

5 

4 



Total incidences = | C | + 2 | D | = 9 . 
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Note that the simplices in A U B form paths: and - a^. 

These paths will be the basis of our algorithms in Chapter IV. 

1.5 Exercises 

5.1 . Show how to construct , for any C C R n that is convex but not compact , a 
continuous function f: C •* C without fixed points. 

5 . 2 [69]. Let g: B n R n be continuous. Show that either 

(i) there is an x* £ B n with g(x*) =0; or 

(ii) there are y*,z* £ 9B n = {x £ R n |jjx|| = 1} with g(y*) = Xy*, g(z*) = yz*, 

X > 0 > y. 

5.3 . Prove the following version of Lebesque's tiling Theorem [42]: For every n, 

there is an e > 0 such that, whenever S n is covered by a finite number of closed 
sets of diameter at most e , there is a point of S n lying in at least n+1 of 
the sets. 

5.4 . For n >_ 2 and e > 0, construct a covering of S n by finitely many closed 
sets of diameter at most e so that no point lies in more than n+1 of the sets. 
(Hint: think of bricks for n = 2 and proceed inductively.) 

3 

5.5 . Show that there is no triangulation of S into more than one regular tetra- 
hedron. (Hint: find the dihedral angle between adjacent faces of a tetrahedron.) 

5.6 . Devise an admissible labelling of the vertices of a triangulation, based on a 
continuous function f, so that the factor n in the inequality of 3.7 can be 
eliminated. 

5.7 . Prove Sperner’s lemma directly from Brouwer’s theorem. (Hint: if each vertex 

of the triangulation is assigned an image in S n and each simplex of the triangula- 
tion is mapped linearly, the result is a continuous function. If a vertex labelled 
i is mapped to v 1 fixed points do not necessarily lie in completely-labelled 
simplices — modify this rule to obtain the desired function.) 

5.8 . Let D n denote the unit cube (x € R ^| x and let f : D n -*■ D n be contin- 

uous . Define the following sets : 



Cq = {x £ D n |f^(x) > x_. or x_. = 0 for all j £ N} and 
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= (x ^ D n |i = min{j|f\(x) <_ x_. > 0 and x_. is maximum satisfying 



this condition}} for each i £ N. 



(i) Show that the C_. , j £ N Q , partition D . Let CL = C_. for each j £ N Q . 
Show that a point in fl. VT C. is a fixed point of f. 

(ii) State and prove a lemma for D analogous to the K-K-M lemma for S . Use 
your lemma to prove the existence of a fixed point of f. (The sets in your lemma 
need not be based on a function f, but the properties they must satisfy can be 



motivated by the C^'s.) (Hint: Construct a homeomorphism h of D onto an 
n-simplex S so that each face of D n is mapped into a subset of a face of S. 



Then apply the K-K-M lemma to S . ) 

(iii) State and prove a result analogous to 3.7 for D n . (This exercise, and its 
continuations in III. 6.4, IV. 7. 5, are related to the "cubical algorithms" of Allgower 
[1] and Jeppsen [31].) 

5.9. Prove the following extension of Brouwer's Theorem: Let T n = 

(x £ R n+1 |v^x = l}. Let f : S n T n be continuous, and let f^(x) _> 0 whenever 
x £ S 1 ?. Then f has a fixed point x* £ S n . 




CHAPTER II: SOME APPLICATIONS OF BROUWER’S THEOREM 



We will return later to the problem of constructing algorithms to locate 
completely-labelled simplices. Here we give some examples of the use of Brouwer’s 
theorem in problems of interest. Of course, the existence result is less important 
to us than the construction of a function whose fixed points solve the given problem. 
We will also see some limitations of Brouwer’s theorem that may suggest extensions. 

If the theorem is applicable, finding a function with the required properties is 
usually straightforward. For any point that is not a solution, one generally has a 
good idea of how it should be modified. While an iterative algorithm based on suc- 
cessive modifications may oscillate and fail to converge, defining a function based 
on this modification and finding a fixed point generally solve the problem without 
difficulty. 

II.l A Model of an Exchange Economy . We will consider a very simplified model; the 
assumptions we make here will be partly relaxed in Chapter VI. A good discussion of 
the economic consequences of our assumptions can be found in Arrow and Hahn [2]. 

It is possible at this stage to include production, but only with the following 
very restrictive assumption. Faced with a set of prices, each producer picks a plan 
to maximize his profit; the required assumption would stipulate that his inputs and 
outputs are then continuous functions of the prices. We will be able to incorporate 
the possibility of production without this assumption in Chapter VI; for now we 
ignore production. 

Suppose there are n+1 commodities, indexed by N^. In the general model each 
commodity is a particular good or service at a particular location and date, and it 
is assumed that all possible futures markets exist. In this case economic agents 
make once-and-for-all decisions of all inputs and outputs for all time, and one must 
assume that they have perfect information. However, it is easier to visualize a 
bartering game played just once. 

There are m economic agents, or consumers. Consumer i has an initial endow- 
ment of the commodities, forming a vector w 1 £ R^ + \ a ^- so ^ as a preference rela- 
tion >\ (x >\ z if consumer i prefers or is indifferent to x compared to z). 
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A price system is a nonzero vector p £ R^ +1 , where (if p_. > 0) Pj/Pj units of 
commodity j are necessary to purchase one unit of commodity i. Given such a price 
system p, the budget set of consumer i is B 1 (p) = {x £ Ip^x-w 1 ) <_ 0}. 
Consumer i will choose for consumption a commodity bundle d 1 (p) maximal with 
respect to in B 1 (p). It seems reasonable that d 1 (p) will exist and plausible 

that it will be unique if p > 0, i.e. , all prices are positive (see exercise 4.1). 

We make the simplifying assumption that d 1 (p) is in fact defined and unique for all 
non-zero p >_ 0. (This may be acceptable if consumer i becomes satiated with 
each commodity.) It is clear that d 1 (Xp) = d 1 (p) for say A > 0. We therefore 
normalize the prices to add to one; in other words, we restrict p to S n . Then 
d 1 : S n R^ +1 is a function, and we assume that d 1 is continuous. 

We further suppose that consumer i spends all his income, i.e., that 

p T d 1 (p) = p^w 1 for all p £ S n . This assumption is fairly innocuous, though it 

does conflict somewhat with the satiation referred to above. Then if d(p) = 

d 1 (p) is the aggregate demand and w = w 1 is the total initial endowment, 

T T 

we obtain p d(p) = p w. Let g(p) = d(p) - w be the excess demand vector. Then 

T n 

p g(p) = 0, i.e., the value of the excess demand is zero, for all p € S . This 

equality is known as Walras' Law, after the economist who first formulated this 
model of an economy. Also, since each d 1 was assumed continuous, g is a con- 
tinuous function from S n to R n+ \ 

A price system would be of little use unless it brought into harmony the dis- 
parate actions of all the consumers. Clearly, in order for all desired trades to 
be feasible, it is necessary and sufficient that g(p) be non-positive. 

Thus it is natural to state 



1.1 Definition , p* £ S D is an equilibrium price vector if g(p*) 0. 

Note that if Walras ' Lav; holds and p* is an equilibrium price vector , then 
g^(p*) <_ 0 for all i £ N Q with equality if p* > 0. Thus all markets are in bal- 
ance except possibly those of free goods, which can be in excess supply. 

We now show that under our assumptions an equilibrium price vector p K always 



exists . 
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1.2 Theorem . Let g: S n + R n+ ^ be continuous and satisfy p^g(p) = 0 for 
all p £ S n . Then there is a p* £ S n with g(p*) <_ 0 . 

Proof . We define a continuous function whose fixed points are equilibrium price 
vectors. To construct such a function, consider the following adjustment scheme for 
price vectors that are not equilibria: increase p_^ if g^(p) is positive (com- 
modity i is in excess demand) and decrease p^ if g^(p) is negative (commodity 
i is in excess supply). This classical "tatonnement” may be globally unstable; 
but we will use it as the basis of our function. Consider the function h taking 
p into p + Xg(p), where X > 0. Unfortunately, h does not take S n into itself 
and must be modified. First define h + : S n -*■ R+ +1 by h + (p) = (h*(p) , . . . ,h*(p) ) T , 
with ht(p) = max{0 ,h^(p) } . (We use this notation for the positive part of a vector 
throughout the present chapter.) For p £ S n , h + (p) is nonnegative, but its coor- 
dinates may not sum to one. Finally, define f: S n -* S n by setting f(p) = 
h + (p)/v T h + (p) . 

T + 

We must show that f is well-defined and continuous. Note that v h (p) >_ 0, 

T T T 

with equality only if h (p) = 0. But then h(p) £ 0 and 0 _> p h(p) = p p + p g(p) 

T T n . . 

= p p by Walras’ Law. Since p p > 0 for all p £ S we have a contradiction, 

T + n 

establishing that v h (p) > 0 for all p £ S . Since g is continuous, so are 
h, h + and v T h + (considered as a function from S n to R). Hence v T h + , being 
continuous and positive on the compact set S n , is bounded away from zero. We 
conclude that f is well-defined and continuous. 

Brouwer’s theorem guarantees a fixed point p* of f. Thus h + (p 5 *) = yp* for 
some y > 0. Let I={i|pV>0}. We then have 

(a) for all i £ I, yp? = ht(p*) = h^(p*) = pV + Xg^(p*); and 

(b) for all j £ I, 0 = yp? = ht(p*) >_h.(p*) = Xg_.(p*). 

If i 6 I, g^(p*) = pp* with p = (y-1) A . Then Walras' Law gives 
0 = p* T g(p*) = l T p*g i (p*) = P*PP* = PP* T P*« But p* T p* >0, so p = 0 and 
g^(p*) = 0 for all i € I. Combining this equality with the conclusion of (b), we 
have g(p*) <_ 0 , as required. Q] 
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1.3 Remarks . 

(a) Although the proof that the fixed points of f are equilibrium price vec- 
tors is a little complicated, the construction of f is straightforward. 

(b) The function f used above is taken from Arrow and Hahn; in contrast, 

Scarf [57,58] uses the function f' with f'(p) = (p + g + (p))/(l + v T g + (p)). The 
fixed points of f' are also equilibria, but for computational purposes f is supe- 
rior. Let p* be an equilibrium price vector with p* positive but very small. 

Let p be close to p* with p^ > p^ and g^(p) < 0. Then f'(p) depends on the 

negative value of g^(p) only through Walras' Law, which requires some g^(p) to 
be positive and hence f _! ( p ) < p^ . If p^ is small, this effect can be minimal. 

As a consequence, p can be very close to f'(p) without g(p) being correspond- 
ingly small. On the other hand, under these circumstances one can expect g^(p) to 
be fairly close to zero, and hence f^(p) will depend directly on g^(p). Since the 
scale of f(p) - p is immaterial, X should be chosen small so that f is more 
sensitive to all coordinates of g. To illustrate these ideas, we solved an example 
of Scarf [57; problem 1] using a modification of the algorithm of Chapter VIII. 

Using f with X = 1 required 120 function evaluations; with X = 10 5 only 95 
were necessary. Using f ' , with or without a scale factor multiplying g + , the 
number required was 138, and the final excess demand vector was much larger. 

(c) We proved theorem 1.2 using Brouwer's theorem. Uzawa [70] showed that the 

implication can be reversed: given f: S n S n , define g by g(p) = f(p) - 

T T 

p(p f(p)/p p). Then g satisfies the hypotheses of 1.2. It is easy to verify that 
if g(p*) <_ 0 then f(p*) = p*. 

(d) In Chapter VI we extend our model of the economy to include production and 
demand correspondences (point-to-set mappings) instead of functions. The limitations 
of the present model will help to motivate extensions of Brouwer's theorem in 
Chapter V. 

II. 2 Unconstrained Optimization 

Just as in section 1 we had to make simplifying assumptions to allow the use of 
Brouwer's theorem, here we will have to ignore constraints and consider only 
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continuously differentiable functions. 

Suppose we are trying to minimize a convex, continuously differentiable function 
0: R n R over R n . For x* to minimize 9 it is necessary and sufficient that 
V0(x*) = 0. Thus a natural function is f with f(x)=x-V0(x). We need not 
worry about step size or convergence — but unfortunately there seems to be no closed 
n-cell that f takes into itself. 

Let C = lev^ 0 = {x £ R n |0(x) < a} be bounded and assume that there is a 
c £ int C. Then 1.2.2 shows that C is a closed n-cell, but one cannot assert that 
f(C) C C. We will have to resort to a device that is useful for proofs but is not 
much help in computation. (Another approach is given in exercise 4.2.) The diffi- 
culty will be removed in Chapter VI. 

Since C is compact, so is f(C); embed C U f(C) in a large simplex S. We 
have f: C -> S and we want to extend f to, say, f: S -*■ S without creating new 
fixed points. Choose r: S -+ C to satisfy: 



(i) 


r 


is 


continuous ; 






(ii) 


on 


C. 


, r 


is the 


identity; 


and 


(iii) 


if 


X 


£ c 


and A 


i — i 

o 

O' 


Ax • 



A function satisfying (i) and (ii) is called a retraction of S onto C; we will 
give two examples of functions r satisfying (i)-(iii) below. 

Now let f: S S be for. f is continuous and hence has a fixed point x*. 
If x* € C, r(x*) = x*, f(x*) = x*, and V0(x*) = 0. We now show that x* g C 
leads to a contradiction. Let z = r(x*). Then x* = z - V0(z). But for suffi- 
ciently small A, z(A) = Ax* t (l-A)z = z - AV0(z) has 0 ( z(A ) ) < 0(z); hence 
z(A ) £ C, contradicting condition (iii) of r. Thus any fixed point of f is a 
fixed point of f and a minimizer of 0 . 

We have proved that a continuous function achieves its minimum on a compact set 
— a result provable by far simpler means. The result is even less impressive when 
we note that finding the minimizer will be difficult even if we have good algorithms 
for finding fixed points because of the retraction r. Even in the simplest case 
evaluation of r requires a line search. Again we are motivated to extend Brouwer’s 
theorem, to encompass both functions on R n and point-to-set mappings (if 9 is not 
continuously differentiable). 
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1 2 

Two choices for r are r and r , defined below: 

(a) let r^Cx) € C satisfy ||r 1 (x) - x|| 2 <_ j|z-x|| 2 for each z £ C; 

(b) let r 2 (x) = Ax + (l-A)c, where A = max(u £ l|yx + (l-u)c € C}. 

II. 3 The Nonlinear Complementarity Problem 

Let g: ■+ R n be continuous. The nonlinear complementarity problem (NLCP) 

•'•T 

associated with g is to find x* >_ 0 with g(x*) >_ 0 and x 5f g(x*) = 0. If g 
is affine (g(x) = Ax + b), this is the well-known linear complementarity problem — 
see Lemke's survey paper [44]. The NLCP was first studied by Cottle [7]. 

One application of the NLCP is to nonlinear programming. Consider the problem: 
minimize(0(z) | g^(z) £0, i = 1,2 , . . . ,m, z £ R^} . Let n = k + m and x = (z ,w) , 
where w £ R m is a vector of Kuhn-Tucker multipliers. When 0 and all the gJs 

are continuously differentiable, finding a Kuhn-Tucker point is equivalent to solving 

T 

the NLCP associated with h, where h(z,w) = (V0(z) + w Vg(z), -g(z)). 

The following result was proved by More [50]: 

3 . 1 Theorem . Suppose there are w £ R^ and p > flwjj^ such that, for all 
x 0 with ||x = p, max ^£ N (x^-wjg^(x) > 0. Then the NLCP associated with g 
has a solution in C^ = {x £ R^ I 1 |x | [^ p). 

Proof . First we define a function h: R^ -> R^ by h(x) = (x - g(x)) + . Then 

x* is a fixed point of h iff it is a solution to the NLCP. 

Since h is continuous, C’ = C Uh(C) is compact; choose t > 0 so that 

P P 

C’ cz C . Note that C is a closed n-cell. Define the retraction r: C C as 
— x x T p 

follows. If x £ Cp , r(x) = x; otherwise, let A = A(x) = 

min{(p-w. )/(x.-w. ) lx. > p} and r(x) = Ax + (l-A)w. Then r is the radial retrac- 
i i l 1 l 

tion onto C^ with center w; clearly A and hence r are continuous functions of 
x. Now let f ( x ) = h(r(x)); then f: C -* C is continuous. 

T T 

Brouwer’s theorem yields a fixed point x* of f. Suppose that x* £ C . Let 
z* = r(x*) and note that z* >_ 0, ||z* = p. Choose i so that (zV-w^)g^(z*) > 

0. Then either z? > w. , in which case x* > z? > w. > 0 and x? > z* > z? - 

li l—ii— l—ii 

g.(z*); or zl < w. and xV < z? <z?-g.(z*)>0. In either case 

l li l—i i i 
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x? i (z? - g^(z*)) + = f^(x*). Hence x* £ , x* = h(x*) and x* solves the NLCP 

associated with g. [[] 

We chose the existence result above because it is relatively strong and also 
simple to prove using Brouwer's theorem. For other results and applications see 
Karamardian [33], Eaves [11], and Fisher and Gould [59]. 

II. 4 Exercises 

4.1 . Let the binary relation >, defined on R^ +1 , satisfy: for all x,y,z £ R+ +1 , 

(i) x > y or y > x; (ii) x > y and y > z imply x > z; and (iii) = 

{x T £ R n+1 |x' > x} and L = {x' £ R n+1 |x >x'} are closed, and U is convex. 

+ 1 — x 1 — x 

Let w £ R^ + ^ be fixed, and let T n denote {p £ S n |p > 0}. For all p £ T n , 
let B(p) = {x £ R^ +1 |p T (x-w) <_ 0} and D(p) = {x £ B(p)|x > y for all y £ B(p)}. 

Show that D(p) is nonempty and convex for all p £ T n . 

Now suppose further that for all x,y £ R^ + \ x ^ y, either x > (x+y)/2 or 

y > (x+y )/2 is false. Show that D(p) is a singleton (d(p)} and that d: T n + 
n+1 . 

R is continuous. 

+ 

4.2. Let 9 : R n R be convex and continuously differentiable, with C = I ev a 0 = 

{x £ R n | 0 ( x ) < a} bounded and c £ int C. Show that there is some A > 0 so that, 
if f(x) = x - AV0 (x) , f takes C into itself. (See also 1.5.9.) 




CHAPTER III: TRIANGULATIONS 



We now present formally the discussion of triangulations, including a proof of 
the properties used in 1.3.4. We also introduce some particular triangulations that 
are used in Chapter IV for the development of algorithms for finding completely- 
labelled simplices. 

Section 1 discusses affine independence and simplices and Section 2 introduces 
triangulations and derives some of their properties. In Section 3 we define some 
particular triangulations and in Section 4 describe their pivot rules. Section 5 
briefly presents Scarf's primitive sets [57,58]. 

III.l Affine Independence and Simplices . Given points x°,...,x- ) in R m , we can 

form the combination x = . X.x 1 , with all X. £ R. Then x is a linear combi- 

^O l ’ l 

nation of the x^s. Also familiar are nonnegative combinations, where each 



x i € V 


and convex combinations. 


where each 


X. £ R 

l ^ + 


and 


i — i 
ii 

•H 

r< 

o 

II 

*l > *H 
t — J 


We say 


x is an 


affine combination of the 


x^'s if 


to h = 


1 , 


but the X . ' s 
i 


are not 



necessarily nonnegative. 



1.1 Definition . If S CR 1 ”, the affine hull of S, denoted aff(S), is 
the set of all finite affine combinations of members of S. S is an affine sub- 
space if aff(S) = S. The dimension of S i 0 is the dimension of the linear sub- 
space parallel to aff(S), i.e. , aff(S) - {s} for any s £ S. (If S is empty, 
dim(S) = -1.) 

All relative concepts of topology are relative to the affine hull. Thus the 
relative interior of S, rel int S, is {x £ S 1 3 e > 0 s.t. B(x,e) fl aff(S) CS}. 
In particular, we mean the relative boundary of a set S whenever we refer to its 
boundary and denote it by 3S. Thus 3S = S fl (aff(S) ^ S). This notation is con- 
sistent with the definition of 3S n in 1.2.1. 

Just as linearly independent points have the property that none is a linear 
combination of the rest , affinely independent points are defined to have the analo- 
gous property. An alternative statement is 

1.2 Definition , y 0 ,...^ in R™ are affinely independent if X^y 1 = 



0 
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and X^ = 0 imply X^ = 0 for all i. One may easily verify the following 

result . 

1 . 3 Lemma . The following are equivalent for y^,...,y^ € R™: 



(a) 


0 j 

y ,-*-,y 


are affinely independent; 








ri 


... l"" 




(b) 


the matrix 


Y = 


o 

) 


... y5_ 


has rank j+1; 


(c) 


the vectors 


i 

z 


i 

= y 


i-1 

- y 


1 < i < j , are linearly independent ; and 


(d) 


dim({y°, . . . 




= j. 


■ □ 




1.4 


Example . v 


0 


n 

»v 


are affinely independent in R n+1 , and their affine 



hull is (x £ R n+1 |v T x = 1}. 

We can now define simplices. 

1.5 Definition . If y^ , . . . ,y^ in R m are affinely independent , then the 

relative interior of their convex hull, i.e., a = {T-J . X.y 1 |x. > 0 for all i, 

^i=0 1 i 

X^ = 1} is a j -simplex . The vertices of a are y°,...,y''. We write a = 
(y° , . . . ). Note that a is relatively open. Denote by a = [y° , . . . ,y^ ] the 
closure of a — a is called a closed j -simplex . 

Note: From 1.3 no simplex in R m has dimension more than m. 

We will not distinguish between a 0-simplex [y] or (y) and its vertex y. 

1.6 Definition . A simplex x is a face of a simplex a if all the vertices 
of i are vertices of a. In particular, a is an improper face of itself. If 
dim t = dim o - 1, x is called a facet of a. If y is the vertex of a which 
is not a vertex of x, x is the facet of a opposite y. 

2 

1.7 Example (m - 2). . y 
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a ~ (y° ,y l ,y 2 ) is the interior of the triangle; t = (y^,y‘*‘) is the open seg- 
ment from y® to y^; i* = (y^) = [y^] is just {y^}. t' and x are faces of 
a , and x T is a facet of x . 

The following important result is easy to prove : 

1 . 8 Lemma . a is partitioned by all the faces of a. ] | 

III. 2 Triangulations . Let C be a convex set in R™. Let A = aff(C) have 
dimension n <_ m. (Note that to triangulate S n we have to consider the case 
n < m. ) 



2.1 Definition . G is a triangulation of C if 

(i) G is a collection of n-simplices; 

(ii) the faces of all the simplices in G partition C; and 

(iii) each x £ C has a neighborhood meeting only a finite number of simplices 
of G. 

Condition (iii) stipulates that G is locally finite. Denote by G^ the col- 
lection of j -simplices that are faces of simplices of G. Thus G n = G. Call mem- 
bers of G^ vertices of G. Denote by G + all faces of simplices of G, i.e., 



U :€N n 



2.2 



(a) 

square : 



(b) 



Examples . 



G = 



Let 



;/ 0 1 2 X , 0 2 3, 

Uy ,y ,y >, \y ,y ,y 




, o 3 4 \ / ° 4 1\, 

<y ,y ,y )> \y ,y ,y » 



triangulates the 



0 < x^ <_ 1}. C can be triangulated as 




follows : 
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This example shows that C need not be closed and, even if C is bounded, G need 
not be a finite collection. 

We summarize the most important properties of triangulations below. Note that 

(a) and (b) give the alternative characterization by closed simplices, and (d) and 
(e) are 1.3.4(a) and (c). 

2 . 3 Theorem . Suppose G is a triangulation of C with dim C = n. Then 

(a) U a 6 G ° = C ' 

(b) If a , a g G and a fl i 0, then fl is the closure of a 

common face t of cr^ and . 



(c) 


If 


DCC is 


compact , D 


meets 


only finitely many simplices 


of 


(d) 


If 


x £ G 11 " 1 , 


then either 










(i) 


T C3C 


and x is a 


facet 


of just one simplex of G; 


or 




(ii! 


) x £ ac 


and x is a 


facet 


of exactly two simplices of 


G. 


(e) 


Let 


D C 3C, 


dim D = n-1, 


and 


D = C f) aff(D). Then G F = 





{ x | t CD, x £ G n 1 ) triangulates D. 

Proof . 

(a) By definition C is the union of all faces of all simplices of G. By 1.8 
a is the union of all faces of a. Thus part (a) follows. 

(b) Let x £ f) a 2 . Then by 1.8 x lies in a unique face of and 

a unique face of Since G is a triangulation, i.e., each 

point of fl lies in a common face of and a^. Let w 0 ,...^ be a list 

of all vertices of such common faces. Then each w 1 is a common vertex of and 

a 2 » Each x £ fl lies in [w°,...,w^] and clearly [w 0 ,...^] C fl 
since the latter is convex. Hence a fl = [w° , . . . jW 3 ] , a common face of 

and o ^ . 

(c) By 2.1(iii), each x £ D has a neighborhood meeting only finitely many 
simplices of G. Clearly, these neighborhoods cover D; hence there is a finite 
subcover. Thus only a finite number of simplices of G meet D. 

(d) Since dim C = n > n-1 = aff(x), there is a point w £ C 'v aff(i). Let 

x be the barycenter of x; i.e., if x = (y^ ,y^ , . . . ,y n ), x = y 1 /n. For 
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l - 1,2,... let x^ = w /2, t (1-1/2, )x. Then x^ £ C for all 2,, x* -* x, and no 

2 , 

x lies in aff(x). There is a neighborhood M of x meeting only finitely many 

i + 

simplices of G but infinitely many x ; hence there is a simplex a in G con- 

2 , 2 , — — 

taining infinitely many x . Since x x, x £ a and x £ x; thus by (b) x 
lies in the closure of a common face of a and x. But x lies in x itself, so 
x is a face of a. Since a / x, a is an n-simplex of G. Thus each x £ G n ^ 
is a facet of at least one a £ G. 



(i) Let x C 3C and assume x is a facet of and a 2 , / cr^. We 

obtain a contradiction. Let = (z^,y^, . . . ,y n ) and = ^w^,y"^, . . . ,y n ). Since 
cr^ is an n-simplex, affCa^) = aff(C) contains w°; hence we have w° = pz^ + 

I i£N X.y 1 with p + l m X. = 1. 

Distinguish three cases: 

Case 1: p > 0 . Then consider c = (l-e)x + ew^ = ew^ + ((l-e)/n)y^ = 

P ez ° + ligx + (l-e)/n)y^. For any 0 < e < 1, c^c^, as shown by the first 

expression. But for sufficiently small e each coefficient in the second expression 



is positive and c £ a^ 9 a contradiction. 

Case 2: p = 0 . Then w^ is an affine combination of y\ . . . ,y n and is 

not an n-simplex. 



Case 3: p < 0 . Every point c in = aff(C) can be expressed uniquely 

Or 1 n 0 

as c = ttw + iKy » y^ = 1. By substituting we find c = pirz + 

^i£N ^i + a l so * If c is close to x, then tt is close to zero and vu 

is close to 1/n, for i £ N. More precisely, x has a neighborhood M in aff(C) 
so that if c 6 M, y. > l/2n and I tt I < l/3nA, where A = max.-.. I A . I . But then 

l — 1 — 101 1 l 1 

any point of M lies in a ^ (if n > 0), x (if tt = 0) or a (if ir < 0), so 
that M C C. This contradicts x € x c 3C. 



(ii) x £ 3C. Then there is an x T £ x with a neighborhood (in aff(C)) con- 

1 2 

tained in C. From this neighborhood we can pick two points w and w with 

12 12 12 
(w +w )/2 = x’ and w ,w £ aff(x). Thus w and w are on opposite sides of 

x. By the construction above we can find two simplices and with x as a 

facet, and it is easy to show that / a^. The arguments of case 1 show that there 

can be no third simplex. 
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(e) G* consists of (n-l)-simplices in G n \ Clearly, all faces of simplices 
of G' are disjoint, and the local finiteness condition holds for G T since it 
holds for G. Thus we need only prove that the simplices of G' and their faces 
cover D. Pick any x £ D; since x £ C, x lies in some p £ G" 1 . Because x 
lies in 3C, it has no neighborhood lying in C so that p cannot be a n-simplex. 
If any vertex of p did not lie in aff(D), then there would have to be vertices 
of p on either side of aff(D). This contradicts D C 3C. Thus each vertex of P 
lies in aff(D) and hence in D; we deduce that p c D. If p is a (n-l)-simplex, 
we are through; otherwise, use the technique of (d) to find a higher-dimensional 
simplex in aff(D) with p as a face. Eventually a (n-l)-simplex in G' is found; 
hence p is a face of a simplex in G’. | | 

III. 3 Some Special Triangulations of R n and S n . For consistency of notation we 
sometimes use names for triangulations that have been applied elsewhere in the 
literature . 

3.1 The Triangulation of R n . This triangulation is often called K or 
"Kuhn's triangulation." Kuhn [37] credits it to Tucker [43, page 140], but it 
appears earlier in Freudenthal [23]. Freudenthal constructed K^ to answer a ques- 
tion of Brouwer: do there exist regular triangulations of a simplex of arbitrary 
mesh whose simplices do not become (in a precise sense) "arbitrarily long and 
skinny"? Whether Brouwer was considering algorithms to compute fixed points is a 
fascinating question. 

Note that our K-^ is not the same as Eaves' in [14]. 
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Let kJ = Z n = {y £ R n | y £ Z for each i £ N}. If y° £ and it is a 
permutation of N, then denote by k^(y^, it) the n-simplex {y^,...,y n ) where y^ = 
y^ ^ + u* (i) for each i £ N. (1.3(c) shows that k^y 0 ,^) is an n-simplex.) 
Finally, let be the collection of all such k^(y^,Tr). Whenever we set 

^(y 0 ,") = <y° . . ,y n >, we suppose the y^s ordered as above, and similarly for 
the other triangulations below. 

3 . 2 Lemma . is a triangulation of R n . 

Proof . Clearly, is a collection of n-simplices and any point of R n has 

a neighborhood meeting only finitely many simplices of K^. We need only prove that 
the simplices of K* partition R n . Let x be an arbitrary point of R n . For each 
i £ N, let y? = |_x^J where, for any real X, L X J denotes the greatest integer 
not exceeding X. Then y^ £ K^. Let z = y^ - x; we have 0 <_ z <_ u. 

Let it be a permutation of N such that 1 >_ Z 7T ( 1 ) . • . >_ z 7T ( n ) Denote 

the terms above by a_= 1, a, , . . . ,a , a . = 0. For i £ N A let 3. = a. - a. . 

J 0 1 n n+1 0 li l+l 



Then each 



3. > 0 and 7. 3. = a - a 

l — i£N^ l 0 i 



n+1 



= 1 . 



Let a = (y°, . . . ,y n ) = k^(y^, it) and consider 3^y 1 . ^ ave 



6 i yl = y0 + B i (yl - y 1 = y + I ieN s i ( ^=i (y] - 



0 

= y + 


^•i£N B i 


0 

= y + 


^j€N (y 


0 


V IT' 


= y + 


^j€N u 



j - y^ 1 )) 






l. u n(j) o. 
^eN D 



tt( j ) 



TT(j) 



= y + z = x. 



Thus x £ a. Using 1.8, we deduce that the simplices of cover R . 

We must now show that all these faces are disjoint. Again let x £ R n be arbi- 
trary. We may assume that x £ a for a = (y^,...,y n ) = k^(y^,7r) £ K^, say x = 

y_T 3.y\ 3. > o for i £ N. and 3. = 1. Then x lies in a face of a 

l i£N q i j 5 l — 0 l i£N q 1 

whose vertices are those y^ for ^£L={££Nq| 3^>0}. We show below how each 

l 0 

y , l £ L, can be generated from x independently of y and tt . Thus these 



same vertices are found for any simplex of whose closure contains x. Hence x 
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lies in just one face of a simplex of and the argument is complete. 

For any z £ a, we have z = 3^(z)y 1 , say, with 3^(z) >_ 0 for i £ N Q 

and 3.(z) = 1. We then form a.(z) = Y. . 3.(z) for i = 0,1,..., n+1 (con- 

l i€N q i l h = i ] 

ventionally ct^ + ^(z) = 0). By the argument above in reverse, we have 1 = <*q(z) 
>...>a .(z)=0 with a.(z) = z ,. N - y ... for i f N. This is true for x 

j i l l 

and also for each y , l £ L. Then we find that for i £ L, (a^(y ^*‘ ,a n+l^ ^ 

has the form (1 ,1 , . . . ,1 ,0 ,0 , . . . ,0 ) with the last 1 in position £. We can now, 

l l 

with our knowledge of the M a-vectors" of x and each y , construct all the y , 
l £ L, directly from x. 



For each y £ (0,1], let y(y) (one of the y ? s) be defined by 



y 1 (y) = 



Lx. J + 1 if x. - Lx. J > y 

l li — 

Lx.J if x. - Lx.J < y. 
i ii 



j i 

It is clear that each y is of this form; also there will be one distinct 
y(y) for each strict inequality in a Q (x) ^ ••• :L a n+ q( x )’ and t ^ lus l L l an tota l- 
Noting that the y(y)'s can be determined from x alone without reference to 
y° or tt, we have completed the proof of the lemma. □ 

3 . 3 Example (n = 4). Let x = (2-^-, 1^-, -2^-, ^-) T . Let us find all simplices 

T 

of whose closures contain x. We have y = (2,1, -3,0) and z = x-y = 

113 1 1 2 

(— , — , — , — ) . We can choose tt = (3, 1,2, 4) or tt = (3, 2, 1,4). In either case 
a = (1, — , y, i, 0) and 3 = (^-, 0, i, ^0 . The two possible simplices are 

<y°,. . . ,y 4 ) = k 1 (y,7r 1 ) and (z°,...,z 4 ) = k^y,^) where y° = z° = (2,1,-3,0) T , 

y 1 = 2 1 = (2,1,-2,0) T , y 2 = (3,1,-2,0) T , z = (2,2,-2,0) T , y 3 = z 3 = (3,2,-2,0) T 
and y - - (3,2,-2,l) T . The reader can confirm that x = \ g.y 1 = I B.Z 1 - Note 

that x lies in the common face (y° ,y\y 3 ,y 4 ). 

We can also proceed directly from x to determine the vertices of the face of 
containing x, as in the last part of the proof. For 0 < y < ^, y(y) = 

(3,2,-2,l) T = yt For Y <i y(y) = (3,2,-2,0) T = y 3 . For y < y <_ |, 

y(y) = (2,1,-2,0) T = y 1 . Finally, for y < y £ 1, y(y) = (2,1, -3,0) = y°. 
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3.4. Note that the crucial part of the proof of 3.2 was obtaining the alterna- 



tive description of k^Cy,^) as those x with 



1 > x 



7r(l) 



Ml) - 



Mn) 



- y. 



ir(n) — 



> 0 . 



(*) 



It is easy to see that k^Cy,^) consists of those x satisfying (*) with 
strict inequalities throughout. By requiring all inequalities but one to be strict 
and replacing the remaining inequality by an equality, we obtain the set of x lying 
in a particular facet of (*). For this reason we refer to (*) as a facetal descrip- 
tion of k^Cyj'rr). From (*) we notice that x £ R n lies in a simplex of K* of 
dimension less than n if and only if some x^ is integer or some (x^ - x_.) is 
integer. 

3.5 The Triangulation of R n . A second triangulation of R n is J^, origi- 
nally due to Tucker [43, page 140] and also found in Whitney [72, page 358]. See 
also [64]. 

Let = Z n . Distinguish certain members of (to be called "central ver- 
tices") by denoting by J^ C the set {y € jJ|y^ f° r eac ^ i € N). If 

y° £ jJ C , it is a permutation of N, and s € R n is a sign vector (each s^ is 
±1), then denote by j 1 (y°,TT,s) the n-simplex (y°,...,y n ) where y 1 = y 1 1 + 

s for each i £ N. Let J. be the collection of all such j,(y°,ir,s). 

ttU) 1 1 

We leave it as an exercise to prove that is a triangulation. The proof 

follows that for K ; the key string of inequalities analogous to (*) is here 



1 ^ Wsu) 



Mi)' 



— S Ti(n)^ X n(n) 



- ^(n)) 



> 0. 



(**) 



3.6 Obtaining triangulations of S n . For any C C R n and 6 £ R, 5C denotes 

{6c|c € C}. For any family G of subsets of R n and any 0 i 5 € let 6G 

denote {6C|C 6 G}; if G is a triangulation of D, 6G is a triangulation of 6D. 

In particular, we now have triangulations 5K^ and 6J^ of R . Now recall that 

mesh G = sup „ diam C. Hence we easily obtain from the definitions that 
p *X€G p 

mesh^SK^ = mesh^SJ^ = Vn" |($| and mesh^SK^ = mesh^S = |5|. Thus we have 
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triangulations of R n of arbitrary mesh. These generate triangulations of S n as 
follows. Denote by C n the set {x £ R n 1 1 > x > . . . ^ x >0} as in 1.2.3. We 
first triangulate C n and then use the method of 1.2.3 to obtain triangulations of 



l \> /V; 

Choose 0 < m £ Z and let (J^) be the collection of n-simplices of 
(J^) meeting mC n = (x € R n |m ^ * n 0}. Let x € mC n lie in a € K 

(a £ J ). As in the proof of 3.2, each vertex of o can be obtained from x by 

sliding each coordinate of x up or down to an adjacent integer. Hence each such 
n 

vertex also lies in mC , and it is easy to deduce that and J triangulate 

n - 1^ - I'v n 

mC , and that m and m triangulate C . Finally, using the affine 

homeomorphism h of 1.2.3, we obtain the triangulations K^Cm) and J 2 (m) of S n . 



3.7 Definition . Let Q denote the (n+l)*n matrix 




and q J 



its 



0 



-1 

+1 



jth column for j € N. Let K 2 (m) = {y € S n |my^ £ Z} for each i £ N Q . If 
y° G K 2 (m) and tt is a permutation of N, let cr = (y^,...,y n ) where y 1 = y 1 1 
t m ^q 77 ^ for each i £ N. If a c s n , we write cr = k 2 (y°,ir) (m is implicit). 
Finally, let ]< 2 (m) be the collection of all such k 2 (y°,Tr). 



3.8 Definition . Let J 2 (m) = K 2 (m) and J^Cm) = (y € J 2 (m)|my^ is even, 

1 i < n and my^ is odd}. If y° € J 2 C (m), 77 is a permutation of N, and 

- _n . . , / 0 n \ , i i-1 -1 77 (i) 

s € R is a sign vector, let cr = \y , . . . ,y ) where y = y +m s^^q 

for each i € N. If a cs 11 , we write a = j 2 (y°,TT,s). Finally, let J 2 (m) be 

the collection of all such j 2 (y°, 7T ,s). 



3 . 9 Lemma . K 2 (m) and J 2 (m) are triangulations of S n . mesh oo K 2 (m) = m 1 , 

mesh oo J 2 (m) = 2m 1 , mesh 2 K 2 (m) = m 1 /n+l if n is odd, m 1 v'n’ if n is even, 
and mesh 2 J 2 (m) = m 1 /4n-2 . 



3.9 and 2.3(d), (e) establish the earlier properties I.3.4(a)-(c) used in the 
proof of Sperner’s lemma. The proof of Brouwer’s theorem is finally complete! 



3 . 10 Remarks . For m a power of 2, J 2 (m) is due to Whitney [72, pp. 358-60]. 
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Note that we have not mentioned the triangulation most familiar to topologists, that 
of iterated barycentric subdivision. The example below shows clearly how this tri- 
angulation progresses. 




The reason for this seemingly glaring omission is that this triangulation exhibits 
the inefficient behavior of creating long, skinny simplices — an algorithm generates 
a large number of simplices ending in one of relatively large diameter, which gives 
a poor approximation of a fixed point. Another disadvantage is that it is difficult 
to obtain the adjacent simplex to a given one when a certain vertex is dropped, and 
even harder to find the new vertex. Shapley has given an algorithm [59] based on 
this triangulation, but there is little to recommend its use. 



3.11 Examples (n = 2). 








y<£.£.°) T ,<l,2)) y<£.H' ,T * (2 > 1)) j 2 ((f,o i) T ,(2,l),(-l,-l ) T ) / 

i 2 (( M’7 )T ’ (2,1) ’ (+1, ' 1)1) 

III. 4 Pivot Rules . All our algorithms will generate a sequence of simplices 



» • • • of a triangulation G with a. and a. adjacent: we obtain a. 

^ 1 1+1 l+l 



VV- 

from a i by dropping a selected vertex y _ of ck and adding a new vertex y + . 
The rules for obtaining from a and y" are called the pivot rules of G. 

Since typically hundreds of these "pivots" are necessary to obtain a sufficiently 
accurate fixed point, it is obvious that the pivot rules should be simple. 



4.1 K^. Let a - <y°,...,y n ) = k^Cy 0 ,^) be given. We wish to obtain x = 

/ 0 n. 0 s i 

\z ,...,z ; = k^(z ,p), with all vertices of c except y vertices of x. The 
table below shows how z° and p depend on y°, xr , and i — from this table it 
is easy to obtain each z and in particular the new vertex of x. 





0 






z 


P 


i = 0 




0 tt(1) 

y + u 


(tt(2) , . . . , tt ( n ) ,tt(1) ) 


n 

•H 

O 


0 

y 


(ir(l) , . . . , tt ( i+ 1 ) , tt ( i ) , . . . ,rr(n) ) 


i = n 


0 7r(n) 

y - u 


( TT(n ) ,tt( 1) , . . . ,Tr(n-l) ) 



tt(1) 



and 



ir(n) 

u 



•n-(l) 



tr(n) 



For 6K 1 replace u 



by 6u 



and 6u 
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4.2 K 2 (m). Let a = k 2 (y°,TT) = (y°,... ,y n ) and x = k 2 (z°,p) 

vertices of a except y 1 . Then and p are obtained from y^ 

according to the table above, except that m and m ^"q 71 ^ 11 ^ 

j ’ n '(n) 

and u respectively. 



contain all 

, tt, and i 

, xr(l) 

replace u 



4.3 J^. Let a = j^(y^,ir,s) = (y°,...,y n ) and x = j (z°,p,t) contain all 
vertices of a except y 1 . Then z^, p and t are obtained from y°, tt, s, and 
i as shown in the table below. 





0 

z 


p 


t 


i = 0 


y + 2s ir(l) U (1> 


TT 


s - 2s , y (1) 

TT(1) 


0 < i < n 


y 


( TT ( 1 ) , . 


. , tt ( i+ l),Tr(i),. 


. . ,TT(n)) 


s 


i = n 


y 


TT 


s - 2s , y (n) 

TT(n) 



For 6J replace u 17 ^ in the top left-hand corner by Su 77 ^. 

4.4 J 2 (m). Let a = j 2 (y°,TT,s) = (y°,...,y n ) and x = j 2 (z°,p,t) contain 
all vertices of a except y 1 . Then z^ , p, and t are obtained from y^, tt, s, 
and i according to the table above, except that m ^q 77 ^^ replaces u 77 ^ 1 ^ in the 
top left-hand corner. 

4 . 5 Remarks . One can check that the simplices above do share the appropriate 
vertices. The pivot rules can be motivated by examining the strings of inequalities 
(*) and (**) of 3.4 and 3.5. Moving between adjacent simplices is equivalent to 
crossing their common facet. Thus two terms in (*) or (**) change position, or an 
extreme term reappears on the other extreme. Clearly, the tables above reflect this 
behavior. 

The pivot rules for iterated barycentric subdivision have been described by 
Shapley [59]. 
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III. 5 Scarf's Primitive Sets . All our algorithms will be based on triangulations. 
However, the first fixed-point algorithm of this type used an alternative notion of 
Scarf [57,58], that of primitive set. We will see in the next chapter that the 
essential property of a triangulation is that of 2 . 3(d) (ii) ; if a is an n-simplex 
of G, t a facet of a, and x £ 9C, then there is precisely one other n-simplex 
o' of G with x as a facet. The notion of primitive set also distinguishes sets 
of n+1 points (like the vertices of a simplex) satisfying a similar condition. 

Let S n denote the set (x £ R n+1 |v T x =1, x <_2v}; S n is a closed n-simplex 
with vertices y^, i £ N Q , where y^ = 2v - (2n+l)v 1 . Pick points y n+1 ,...,y^ 
in S n such that for all i £ N Q and l > n, m > n with l i m, y* i y™. 



Z Ok 

5.1 Definition. A set of n+1 points (y |il £ L} from (y , . . . ,y } form 



a primitive set if there is no m, n < m <_ k, with y^ > min^^ Y^ f or each 
i £N 0 . 



5 . 2 Example (n = 2). We have y° = (-3,2,2) T , y 1 = (2,-3,2) T and 
y 2 = (2 ,2 ,-3) T . Let y 3 = (.9,0,.1) T , y 4 = (.1,.9,0) T , y 5 = (0,.1,.9) T and 

y 6 = ( .05,.5,.45) T , y 7 = ( . 45 , . 05 , . 5 ) T , y 8 = ( . 5 , . 45 , . 05 ) T . The primitive sets 

r 1 2 3, , 0 4 6, , 0 5 6, ,37 8,^ , 

are {y ,y ,y } etc., ly ,y ,y }, ly ,y ,y } etc., iy ,y ,y } etc. and 

{y 8 ,y 7 ,y 8 }. (Here "etc." means other triples obtained by the permutation 

, 0 1 2 W 3 4 5 W 6 7 8. . 

(y y y )(y y y )(y y y ) . ) 

One can see the close relation of these primitive sets to the triangulation 
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5 . 3 Lemma (Scarf [57]). 

(a) Let {y^|& £ L) be a primitive set and l £ L. If {y^ | Jl £ L ^ }} £ 

(y 1 1 i £ N q }, then there is a unique J l + with {y^ | Jl £(L^{A})U{S- + }} a 
primitive set. 

(b) If (y m |m £ M} c (y^i £ N Q ) and |m| = n, then there is a unique m + with 

{y m |m ( M U {m + }} a primitive set. [ | 

5 . 4 Remarks . In the next chapter we briefly describe how Scarf’s algorithm 
uses primitive sets. Let us now compare primitive sets and triangulations. First, 
triangulations are essential for the more sophisticated algorithms we discuss in 
Chapter IX. However, primitive sets have the apparent advantage that the set of 
points y n+ \ . . . ,yk determines the primitive sets, while there are many triangula- 
tions having these as vertices. Thus if a particular arrangement of vertices seems 
desirable, much less memory is required for primitive sets rather than triangulations. 
On the other hand, the replacement operation of 5.3(a) necessitates a search through 
y n+ \ . . . ,yk unless these points are chosen with sufficient regularity. The almost 
universally used method of choosing these points leads to a direct correspondence 
between primitive sets and simplices of K^(m) for some m > 0. The method is due 
to Hansen; see Chapter 7 of Scarf [58]. 

III. 6 Exercises . 

6.1. Prove 1.3. 

6.2. Let C be as in Section 2. Let G satisfy (i) and (iii) of 2.1, and (a) and 
(b) of 2.3. Show that G triangulates C. 

6.3. Prove that ^ triangulates R n . Show that x £ R n lies in a simplex of J* 
of dimension less than n if and only if some x^ is an even integer, some x^ - x_. 
is an even integer, or some x^ + x_. is an even integer. 

6.4. Prove that the simplices of nT^K^ (or of m 1 J^) lying in D n = 

(x £ R n |0 <_ x <_ u} triangulates D n , for any integer m > 0. 

6 .5. Prove that ^(m) triangulates S n . 

6 . 6. Prove 5.3. 

6 .7. Let (x° jy 1 , . . . ,y n ) and (z° ,y^ , . . . ,y n ) be non-intersecting n-simplices lying 
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in an affine subspace of dimension n. Let X = 

“l 1 ... 1 “ 



Z = 



0 1 
z y 



0 1 

x y 



and 



Show that the determinants of X and Z have opposite 



signs . 




CHAPTER IV: ALGORITHMS TO FIND COMPLETELY-LABELLED SIMPLICES 



While the proof of Sperner's lemma given in 1.4 does not suggest any method to 
find completely-labelled simplices, we noted in example 1.4.3 that n-simplices of G 
with the labels 0,1,..., n-1 formed paths. Cohen [6] gave a proof of Sperner's 
lemma based on these paths; we present his argument in Section 1. Cohen's proof is 
inductive — we still have the problem of how to start. Two possible methods will be 
given in Section 2. Sections 3 and 4 show how these methods can be implemented 
using K 2 (m). Section 5 describes Scarf's algorithm, and Section 6 illustrates the 
algorithms with an example. 



IV. 1 The Graph T . Suppose we have a triangulation G of S with the vertices 
of G admissibly labelled, i.e., no vertex in S^ has the label i. We call mem- 
bers of G and G n ^ merely n-simplices and (n-l)-simplices . An n-simplex is 
completely-labelled (c.l.) if its vertices carry all labels in N Q , while an n- 
or (n-l)-simplex is almost completely labelled (a. c.l.) if its vertices carry all 
the labels 0,1,..., n-1. 

We can now form the following graph. 



1.1 Definition. Let the nodes of T be a. c.l. n-simplices and a. c.l. (n-1)- 
n 

simplices lying in S^. Two nodes of are adjacent if one is a face of the 

other or if they share an a. c.l. face. 



2 

v 
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The nodes of T are the heavy dots and its edges are the double lines. (The 
dot is placed on the barycenter of the corresponding simplex.) 

We now analyze the degree of each node of , i.e., the number of adjacent 

nodes . 

Case 1 . t is an a.c.l. (n-l)-simplex in S^. By III.2.3d(i), t is a face 
of just one n-simplex a and clearly a is a.c.l. Hence x has degree one. (See 
x in 1. 2 . ) 

Case 2 . a is a c.l. n-simplex. Then a has just one a.c.l. facet x, that 

opposite the vertex of a labelled n. If x C S^, x is a facet only of a 

(III.2.3d(i)) ; a is adjacent to x and to no other node of T^. If x £ S^, 
then since the labelling is admissible we have x £ 3S n and x is a facet of 
exactly one other n-simplex o' by III. 2. 3d(ii) . Clearly o' is a.c.l. Thus in 
either case a has degree one. (See c and a ^ in 1.2.) 

Case 3 . o is an a.c.l. but not c.l. n-simplex. Since the n+1 vertices of 

1 2 

o have n labels, there is precisely one pair of vertices, y and y , with 
the same label. Thus a has precisely two a.c.l. faces x and x^ opposite y 1 

and y^ , respectively. As in case 2, each x^ either lies in S^ and is a node 
of T or leads to another a.c.l. n-simplex ck. In any case a has degree two. 

(See a and in 1.2.) 

By combining these cases we immediately obtain 

1 . 3 Theorem . Each connected component of has one of the following forms: 

(i) A simple circuit whose nodes are a.c.l. but not c.l. n-simplices. 

(ii) A simple path whose intermediate nodes are a.c.l. but not c.l. n-simplices, 
and each of whose endpoints is either 

(a) a c.l. n-simplex or 

(b) an a.c.l. (n-l)-simplex in S^. | | 

Example 1.2 shows that all four forms can arise. 

Since the number of endpoints of a path is two, the total number of c.l. 
n-simplices and a.c.l. (n-l)-simplices in S^ is even. Clearly this proves the 
inductive step of Sperner’s lemma. 
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IV. 2 How can we find a c.l. n-simplex using First note that even if we have an 

a.c.l. (n-1) -simplex in S^J, tracing the path in from this endpoint may take us 

back to S^, as in the leftmost path in 1.2. The only way to be sure of avoiding 
this behavior is to have only one a.c.l. (n-l)-simplex in S^. 

The first method, to be discussed in Section 3, ensures that this is true by 
adding an artificial layer to S^. The second method, described in Section 4, 
mimics the inductive proof by linking together all the graphs I\ , j £ N. Before 
we give details of these methods using K^Cm), we can give an intuitive idea of 
their operation. 




Note that there is now a path from the unique a.c.l. 1-simplex to a c.l. 2-simplex. 
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O 1 

Now associate with S x and construct 

1 
v 



1 
V 

» c ln 0 

— v — to a c.l. 

2-simplex. 

IV. 3 The Artificial Start Algorithm (Kuhn [38]). We have a given function 

f: S n ■+ S n . We obtain an approximate fixed point of f by finding a completely 

labelled simplex of K^Cm) using the standard labelling. To do this we first show 

how to extend K (m) to a triangulation of S n with an extra '’layer" below S^. 

The arguments of III. 3.6 show that the simplices of that meet mC n = 

{x £ R n |m >_ x >_ -1} triangulate mCp; hence those of m 1 K^ that meet 

C n = (x£R n |l > Xl > ... >_ x r >_ -m triangulate C n . By applying the homeomor- 

phism h, we obtain a triangulation K 2 (m) of S n = {x £ R n+1 |v T x =1, x^ 0, 

i = 0,...,n-l, x > -m" 1 }. S n is S n with an extra layer added to S n , as 
n — n 

described informally in Section 2. 

Each vertex of K 2 (m) lies in S n or has its last coordinate -m . If 
y £ K^Cm) lies in S n , we label it min{j £ [ f ^ ( y ) y_. > 0}; if y £ S , we 

label it min{j € N Q | y ^ = max^h The latter labelling roughly corresponds to 
extending the function f to S n so that f(x) = (n+1) ^v for 



Then the union of T and T 2 forms a graph T : 
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X £ S n = {x e S n |x = -m 1 }. 
n 1 n 

If we let S? = {x £ S n |x i = 0} for i = 0,...,n-l, we clearly have an admis- 

sible labelling, in the sense that no vertex in has the label i. Also, it is 

n 

apparent that if a 6 K 2 (m) does not lie in S , then each of its vertices has 
n"^ coordinate 0 or -m ^ and cannot be labelled n. Hence all completely labelled 
simplices of K 2 (m) in fact lie in S n and give approximate fixed points of f. 

We will therefore have the makings of an algorithm if we can establish that there 

^n 

is just one a.c.l. (n-l)-simplex x* in S^. 

Let m = kn for some integer k > 0. Then we have 



3 . 1 Lemma . There is just one a.c.l. (n-l)-simplex x* of ^(m) in 3^. 

Proof . Note that the simplices of ^(m) are defined as in 3.7, except that 
each vertex must lie in S n but not necessarily in S n . We first exhibit an a.c.l 
(n-l)-simplex x* in and then show that it is unique. 



Let y* = ((k+l)m \km \ ...,km \-m and n* - (1 ,2 , . . . ,n-l) . The (n-1)- 

s implex x* = k (y*,Tr*) is defined in the natural way. Its vertices are 



0 

y 


n-1 

• *»y 


with 


0 

y = y~ 


and y 1 


i-1 -1 7r*(i) 

= y + m q 


for 


i = 1,2,... 


,n-l. 


Thus 


i 

y = 


• — 1 

1 

g 


. . ,km 


\(k+l)m 


- 1 w" 1 

,km , . . 


-1 -1 T 

. ,km ,-m ) , with 


the 


(k+l)m ' L 


in the 


.th 

1 


coordinate . 


Then 


y 1 has 


label i 


for each i. Also, 


X* 


is a facet 


of 





K 2 (y*,(l,2, . . . ,n)) . Hence x* is an a.c.l. (n-l)-simplex in sP. 

Now let x = (y ,...,y ) = k 2 (y ,ir) be an a.c.l. (n-l)-simplex in S^. We 

must show x = x*. Let y^ = (a._m \...,a. ,m \-m for each i. Then 

J lO i,n-l 

yri a.. = m+1 = kn+1 for each i. Hence y 1 can have label j only if 
L ]=0 13 

a. . > k+1. 



13 - 



If a Q0 < k+1, each a 
there is some j with a Q _. 
label j. Thus a^ = k+1. 
while if any a Q _. < k then 
In conclusion, y^ must be 
If xr ( 1 ) f 1, then y 1 
be 2, etc., so that it = tt* 



. ^ < k+1 and no vertex has label 0. If a > k+1, 
lO 00 

< k, so that a. . < k for each i, and no vertex has 
i] “ 

A similar argument shows that no a^_. can exceed k+1, 

a > k > a. . for all i and no vertex has label j . 
lO - - 13 

y*. 

has label 0 — a contradiction. Similarly tt ( 2 ) must 
and x = x*. Thus x* is the only a.c.l. simplex 
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of 1 (m) lying in | | 

3.2 Algorithm . We are given f: S n -► S n . Pick m = kn for some integer 

k > 0. Triangulate S with K^Crn), and label the vertices of ]< 2 (m) as above. 

Let x* be the (unique) a.c.l. (n-1 )-simplex of 1< 2 ''“(m) lying in S n . 

o> 

Step 0 : Let be the unique n-simplex of K^Cm) that has x* as a facet. 

Let y + be the vertex of that is not a vertex of x*. Set £ - 1. 

Step 1 : Calculate the label of y + . If it is n, STOP; cr^ is completely 

labelled and yields an approximate fixed point of f. Otherwise, the label of y + 

duplicates that of exactly one other vertex of a , say y . 

Step 2 : Find the simplex a £ +1 that has as vertices all the vertices of a 

except y . Let y + be the vertex of a n . that is not a vertex of a„. Set 

r J J l+l £ 

£ +■ £+1 and return to Step 1. 



3 . 3 Proof of Convergence . First we show that are distinct. Assume 

the contrary, and let a. = a , j < £ be the first duplication. If j > 1, then 
3 £ 

a. is adjacent to a. _,a. . and cr n _ in the graph r . Since no node of r 
3 J 3-1 3+1 £-1 6 ^ n n 

has degree greater than two, two of these must be equal. We cannot have ^ = 

a_. + ^ or a j ]_ = a £ since these would be previous duplications. Thus a j + q = 

a„ .. and j +1 = £-1, i.e., £ = j+2 and a. = a. Now consider the vertex y 

£-1 J » j j -j 3+2 J 

of a . . that is not a vertex of a.. Since the other vertex with the same label 
3+1 3 

as y is dropped to obtain a j + 2 s we ^ ave that y is a vertex of °j + 2 * This is 

a contradiction. On the other hand, if j = 1, then is adjacent to a j + q and 

o. . But since a. has degree one in Y , we have a. .. = a . , yielding the 
£-1 1 6 n’ 3+1 £-1 J 6 

same contradiction. 



Since the simplices generated are distinct and there is a finite number of 

f\f 

simplices of K 2 (m), the algorithm must terminate. If it terminates with in 

Step 2 because does not exist, then has an a.c.l. facet in 3S n . Since 

the labelling is admissible, this facet must lie in S^ and hence be x*. Thus 
o ^ = a , contradicting the fact that all simplices are distinct. 

Hence the algorithm terminates in Step 1 with a completely labelled simplex. Q 
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This is the only time we give a proof of convergence for an algorithm; the 
proof is similar for any algorithm to be described. Usually we merely describe a 
graph with the same kind of property as r . The algorithm is then a formal state- 
ment of: "Start at a uniquely-defined simplex and trace a path in the graph until 

a desired simplex is found." 

3 . 4 Remark . Lemke and Howson [45] first gave a proof of convergence based on 
the novel arguments above in their seminal paper on equilibrium points of bimatrix 
games. The uniqueness of the argument stems from its reliance on purely combina- 
torial reasoning to assure finite convergence — there is no monotonic evolution. 

Such proofs are now standard in complementary pivot theory. Lemke [44] has a 1970 
survey on complementary pivot theory, while Gould and Tolle [26] and the author [63] 
discuss its general application from an abstract viewpoint. 



IV. 4 The Variable Dimension Algorithm . The version of this algorithm we will pre- 
sent is that of Kuhn [39]. Shapley constructed an algorithm independently [59], 
using barycentric subdivision. However, both algorithms are based on the pioneering 
work of Scarf [57] using primitive sets. 



4.1 Definitions . Let = {x £ S | x j +1 = ••• = x n = 0} for 

j = 0,l,...,n-l. s ^j + q) a c i° se d j-simplex naturally corresponding to S' 1 . For 
j £ N, let I\ be the graph whose nodes are j-simplices in 0:r 

simplices in S^_. ^ whose vertices carry the labels 0,1,..., j-1. Two nodes are 
adjacent if one is a face of the other or if they have a common face whose vertices 
carry the labels 0,1,..., j-1. 

For each j £ N, I\ satisfies a result similar to 1.3, where j replaces n, 
S* } replaces S^, and "a.c.l." ("c.l.") is interpreted to mean "whose vertices 
carry all the labels 0,1,..., j-1 ( 0 ,1, . . . , j ) . " 

Let T be the union of , for all j £ N, i.e., the set of nodes (edges) 

of T is the union of the sets of nodes (edges) of all I\ . Then we have 

4.2 Theorem . Each connected component of T is either a simple circuit or a 



simple path from v to a c.l. n- simplex. 
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The proof follows from 1.3. [^J 



A natural triangulation to use is K 2 (m). We need to be able to describe sim- 
plices of K^Cm) in S^ +1 ^. Let f be a permutation of {1,2,. ..,j}. Then, for 
€ K^Cni), denotes the j-simplex (y , . . . ,y^ ) with y* = y* + m q 

for 1 < i <_ j, if all y 1 lie in S (j +1 )- The unique ( j+l)-simplex of K^ +1 (m) 
lying in S^ +2 ^ with k^Cy 0 ,^) as a facet is k 2 (y°,TT ), where 

TT* = (-Tr(l) , . .. ,7T(j),j+l). 

We then have the following. 



4.3 Algorithm . We are given a function f: S R + S D . Triangulate S n with 

K^Cm) and label vertex y with min{j|f^(y) <_ y^ > 0}. 

Step 0 : Let j = 1, y° = v°, a = k 2 (y°,7r) with tt = (1) the permutation 

of {l}. Let y + = y^ = v^ + m ^q\ Set i = 1. 

Step 1 : Calculate the label of y + . If it is j, go to Step 3. Otherwise, 

it duplicates the label of some vertex y of cr^. 

Step 2 : Let t be the facet of opposite y . If t C S^p go to 

Step 4. Otherwise, let he th e unique j-simplex in s ^j +1 ) sharing the facet 

t with a . Let y + be the new vertex of a . Set l £+1 and return to 
i J ^+1 

Step 1. 

Step 3: (Increasing the dimension) We have the j-simplex o ^ = k^(y ,Tr), say, 

with tt a permutation of {l,2,...,j}. The vertices of have all the labels 

0,1,..., j. If j = n, STOP; a ^ is completely labelled and yields an approximate 
fixed point of f. Otherwise, let °£ +1 = k 2 (y°,TT ), where 

tt’ = ( tt ( 1 ) , . . . , tt ( j ) , j+1) , and let y + be the new vertex of u £+1 - Set l **- l+l , 
j «- j+1 and return to Step 1. 

Step 4: (Decreasing the dimension) We have the j-simplex a ^ = k 2 (y ,tt) = 

(y^,...,y^) with tt a permutation of {l,2,...,j}. has a facet t in ^(j)’ 

and it is clear that t must be (y^,...,y^ ) and ^(j) = j- The vertices of 

i have all the labels 0,1,..., j-1. Let a £+1 = k 2 (y°,Tr') with 
tt* = (Tr(l),...,Tr(j-l)), and let y be the vertex of 
l Jl+1, j ^ j-1 and return to Step 2. 



with label j-1. Set 
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This algorithm produces a sequence of simplices of varying dimensions following 
the path of r from to a c.l. n-simplex. 

Alternatively, we can "fill out" each simplex of dimension less than n by 
adding artificial vertices. We then obtain a triangulation of S = 

{x £ R n+1 |v T x =1, x <_ 2v}. (See III. 5.) The new vertex y 1 = 2v - (2n+l)v 1 is 
labelled i-1 (mod n+1). The algorithm is then equivalent to that of Section 3. 
Using the example of Section 2, we obtain the extended triangulation of S n shown 
below. The natural extension of is also shown, and its close relationship with 

F = *1 U r 2 is clear. 



0 2 




IV. 5 Scarf 's Algorithm [57]. The picture above, together with III. 5, suggests how 

i 

an algorithm using primitive sets can be constructed. Let S , y , i £ be as 

i . Z 

in the last section. Label y i-1 mod n+1, for i £ N Q , and label y 

i Z Z 

min{j|f_.(y ) <_ y^. > 0} if Z > n. Start with the unique primitive set including 
{y 1 | i £ N} given by 111.5.3(b). Proceed as in 3.2, using the replacement step of 



111.5.3(a) in each Step 2. 
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Scarf’s algorithm used a different labelling scheme: label y i for each 

i € and label y min{j|fXy ) >_ y^} for l > n. Again a set of close points 
carrying all labels gives an approximate fixed point, and the algorithm is the same. 
Vertgeim [71] discusses some general labelling rules, including both those used here. 



IV. 6 Examples . We take n = 2 so that the progress of iterations can be seen on a 
diagram. We execute the algorithms using the formal descriptions to show how rapidly 
and conveniently the operations can be performed using K^. 



Define 



f: S 2 + S 2 by f(x) 



~.2 .1 .3' 

.3 .4 .3 

__ . 5 .5 . 4 _ 



We will use K^(4) to find 



an approximate fixed point of f. In fact, f has just one fixed point, x* = 
(7 ,11 ,15 ) T /33 . 



6.1 The Artificial Start Algorithm Illustrated . We have m=4=kn=2.2. 
The course of the algorithm is shown in the table below. The operations are per- 
formed from left to right and top to bottom. We omit the denominator of 4 in every 
vertex. 

The picture below shows the simplices generated. 



2 

v 
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Label Table III. 4.1 




a 



9 



is a completely labelled simplex. 



Note that in fact x* lies in a . 

O 



6.2 The Variable-Dimension Algorithm Illustrated . We use K^(4). The table 



shows the course of iterations, as does the diagram following. 
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IV. 7 Exercises 

7.1. Apply the algorithms of Sections 3 and 4 to the labelled triangulation in 
Section 2. 

7.2 . Devise an artificial-start algorithm using a triangulation based on J 2 (m), 
where m = kn and k is a positive even integer. Apply your algorithm to the 
example of Section 6. 

7.3. Devise a variable-dimension algorithm using J^Cni). Apply it to the example 
of Section 6. 

7.4. We have not shown exactly how the triangulation at the end of Section 4 was 
obtained. Let S n = {x £ R n+1 |v T x =1, x <_ 2v} and y 1 = 2v - (2n+l)v 1 for 

i 6 N . If I,J is a partition of N Q with 1^0, let be the n-simplex 

with vertices v 1 , i £ I, and y^ , j € J. Let H be the set of all such a^. 
Prove that H is a triangulation. (Hint: do not use the methods of Chapter III. 
Consider the natural triangulation of the octahedron in 

R n+1 : {x 6 R n+1 |I |x. | = 1}.) Next show how any triangulation G of S n can be 
N 0 1 

'Xj r iin 'v . . 

extended to a triangulation G of S . Prove that G is a tn angulation. 

7.5 . Refer to 1.5.8, III. 6. 4. Triangulate D n with m" 1 ^ or m 1 J 1 * Label a 

vertex y i if y £ c! . Construct an algorithm to give a completely-labelled 
simplex. (Hint: find a set of n labels so that, independent of f, there is a 

unique (n-l)-simplex in the boundary of D n with these labels.) 



CHAPTER V: EXTENSIONS OF BROUWER'S THEOREM 



Chapter II makes clear the need for more general fixed-point theorems and algo- 
rithms to approximate the resulting fixed points. Two extensions are necessary: 
one to relax the continuity of the function involved, and one to relax the require- 
ments on the domain and range of the function. Section 1 discusses upper semi- 
continuity of point-to-set mappings. Section 2 introduces the fundamental concept 
of piecewise-linear approximation and with this tool proves Kakutani’s theorem on 
the existence of fixed points for upper semi-continuous mappings. Section 3 extends 
Kakutani’s theorem in ways that will be useful in applications. 

V.l Upper Semi-Continuity . Let us motivate our extension by considering a simple 

2 

one-dimensional optimization problem. If we want to minimize 0^(x) = x - x > 
we can form f^(x) = x - VQ^x) = x - (2x - 1) = 1 - x. Then [0,1] -► [0,1] 

is continuous and a fixed point gives the minimizer x* = ■— of 0^. Now let 

9 (x) = ]H X " \\ • V0 2 (x) is - i for x < y, for x > j and undefined for 

x = j. Thus f ^ ( x ) = x - V0 2 (x) is j + x f or x < ^ and x ~ \ for x > ^* 

f 2 has no fixed points. Clearly, we must define f 2 (“0 somehow, and thus extend 
the notion of differentiability. 

1.1 Definition . Let 0 : R n -* R be convex. Then h ^ R n is a subgradient 

of 0 at x if h T (z-x) <_ 0(z) - 0(x) for all z £ R n . The subdifferential of 
0 at x, denoted D0(x), is the set of all subgradients of 0 at x. 

It follows from separation theorems that D0(x) is nonempty for all x when 
0 is convex, and clearly D0(x) is convex. If 0 is differentiable at x, 

D0(x) = (V0(x)} , and 0 £ D0(x) if and only if x minimizes 0. The subdiffer- 
ential is therefore a natural extension of the gradient. 

Since D0 is set-valued, we are led to consider point-to-set mappings. For 
our example, we have D0 2 (x) = {- -j} if x < — , = [- 7 ^ + -|- ] and 

D0 2 (x) = {+ y) if x > If we define F(x) = (x) - D9 2 (x), we have F shown by 
the diagram below: 
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2 



We naturally call y a fixed point of F, meaning £ F(y). 

Clearly, we need some conditions on F to guarantee fixed points. Some sort 
of continuity is required, and motivated by the example above we see that we can 
allow F(x) to suddenly increase but not suddenly decrease. 

1.2 Definition . Denote by P(R m ) the set of all subsets of R™. Let 

C C R m and F: C -► P(R P ) be a point-to-set mapping. Then we say F is upper 
semi-continuous (u.s.c.) if 

(i) for all x £ C, F(x) is compact; and 

(ii) for all x £ C, for all e > 0, there is a 6 > 0 such that if 
z £ B(x,6) fl C, F(z) C B(F(x),e). 

If f is a function from C to R P , then F: C -> P(R P ) with F(x) = {f(x)} 
is a point-to-set mapping — we write F = {f}. Note that f is continuous iff {f} 
is u.s.c. 

The most important property of u.s.c. mappings (often taken as a weaker defini- 
tion of upper semi-continuity) is the following: 

1.3 Lemma. If F: C + P(R P ) is u.s.c., {x^} is a sequence of points in C 

k k k 

tending to x*, and {y } is a sequence of points tending to y* with y £ F(x ) 

for each k, then y* £ F(x*)- 

Proof. For any z > 0, we can find k such that ||y* - y e/2 and 

F(x k ) C B(F(x*) , e/2). Thus y* £ B(F(x*), e ) for every z > 0. Since F(x*) is 
compact, y* £ F(x*). Q 
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Clearly, that F is u.s.c. from C to P(C) for some compact convex C C R n 
does not suffice for F to have a fixed point. Indeed, the map that takes every 
point to the empty set is u.s.c. Also, if F(y) in the example above 1.2 is changed 
from [0,1] to {0,1}, then F remains u.s.c. but has no fixed point. These 
examples motivate the restriction that the values of F must be nonempty convex 
sets. If F is not convex-valued, then it can be "convexified" . This is one part 
of the following theorem, which collects all the results we need about composing and 
combining u.s.c. maps. 

1.4 Theorem . In all parts below, unless otherwise specified, all u.s.c. 
mappings are from C C R m to P(R P ). 

(a) If F is u.s.c. and D CC is compact, so is F(D). 

(b) Let F be u.s.c. and G: D C R P -* P(R A ) be u.s.c. with F(C) CD. Then 

H = GF : C P(R £ ) is u.s.c., where H(x) = G(F(x)). 

(c) Let F^ be u.s.c., i = l,2,...,k and F: C -*■ P(R^) be defined by 

F(x) = if , F. (x). Then F is u.s.c. 
i=l l 

(d) Let F^ be u.s.c., i = l,2,...,k, and F: C P(R^) be defined by 

F(x) = F^(x) + ... + F^(x). Then F is u.s.c. 

(e) Let F be u.s.c. and C be closed with C CD. Then G: D -* P(R^) is 

u.s.c., where G(x) = F(x) if x £ C, <j> if x £ D ^ C. 

(f) If F is u.s.c., then conv F: C + P(R^) is u.s.c., where (conv F)(x) = 

conv(F(x) ) . 

p. 

(g) If F^:CCR n ->P(R 1 ) is u.s.c., i = l,2,...,k, and p = £ p^, , then 
F is u.s.c., where F(x) = F^(x) x ... x F^(x). 

Proof . For (a), (b), (c) and (g), see Theorems 3, 1', 3* and 4' (pages 116-120) 

12 k 

of Berge [3]. Part (d) follows from (g) and (b), with G(z ,z ,...,z ) = 

1 k 

z + ... + z . Part (e) is trivial from 1.2. We now prove (f). 

Given x £ C and e > 0, determine 6 > 0 so that z £ B(x,<5) 0 C implies 

F(z) C B(F(x) ,e ) ) . Let g £ conv F(z), so that g = A^g 1 , >_ 0 and 

A = 1 for some g^£F(z), i=l,...,k. For each , there is an 

f 1 £ F(x) with Hg 1 - f 1 1| <_ e. Since the Euclidean norm is a convex function, 
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rk i 

||g - f|| 2 1 where f = A^f £ conv F(x). Hence conv F(z) c B(conv F(x),e). 

Clearly, conv F(x) is compact for any x £ C, and thus conv F is u.s.c. 

We generally use 1.4 as follows. In different (closed) regions, we have dif- 
ferent u.s.c. maps (often functions) corresponding to a modification of the given 
point. These maps are extended to R m by 1.4(e) and joined together by 1.4(c), 
and the result M convexified M by (f). 

m* 

It is convenient to denote by R the set of all nonempty convex subsets of 



V.2 Kakutani's Theorem and Piecewise Linear Approximations . In 1941 Kakutani [32] 
generalized Brouwer's theorem to u.s.c. maps. He proved 

2.1 Theorem . Let C be a compact convex subset of R m , and F: C -► C* be 
u.s.c. Then F has a fixed point x*, i.e., x* £ F(x*). 

To prove 2.1 we use the fundamental notion of a piecewise linear approximation 
to F. 



2.2 Definition . Let OCR be convex, with dim C = n, and let G be a 
triangulation of C. Let F: C + P(R^) be a point-to-set mapping whose values are 
nonempty. For each vertex y of G, pick f(y) £ F(y). Each x £ C lies in a 
unique simplex of G + . Thus we can write x = A^y 1 , A^ >_ 0, A^ = 1, 

where (y^,...,y n ) £ G and the pairs A^, y 1 for A^ > 0 are unique. Then set 
f(x) = ]j\^ A^f(y 1 ). The function f: C + R^ is thus well-defined; we call f a 

piecewise linear (p.l.) approximation to F with respect to G. 



2 . 3 Lemma . A piecewise linear approximation is a continuous function. 

Proof . An easy exercise using the local finiteness of G. | | 

Proof of 2.1 . We will follow Eaves' proof [12], which is closely based on 
Kakutani's original proof. If dim C = n < m, we can map aff(C) homeomorphically 
onto R n . We can therefore assume without loss of generality that m = n. Pick 
c £ int C. Since C is compact, we can embed it in a closed n-simplex S. Extend 
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F to S as follows. Let F^x) = F(x) if x £ C, F^x) = 0 if x £ S ^ C. Let 
F 2 (x) = {c} if x £ S ^ int C, 0 if x £ int C. Finally, let F Q (x) = 
conv(F^(x) U F^Cx)). By 1.4(e), (c) and (f), F^ : S *> S* is u.s.c. 

Further, if x* is a fixed point of F Q , then certainly x* £ C (otherwise, 

x* £ Fq(x*) = {c} cC, a contradiction). If x* £ int C, then x* £ F^(x*) = 

F^(x*) = F(x*). So assume x* £ 9C. Then there is f £ F^(x*) = F(x*) and 
X £ [0,1] with x* = Xc + (l-X)f, since F(x*) is convex. If X > 0, then 

since c £ int C and f £ C, x* £ int C; thus X = 0, x* = f, and x* £ F(x*). 

Hence x* is a fixed point of F. 

We now prove that F Q has a fixed point. Let G^, k = 1,2,..., be a sequence 
of triangulations of S with mesh -► 0. (Such triangulations exist; for example, 
K 2 (m) can be mapped by a function taking S n onto S barycentrically . ) For each 
k, let f^ be a p.l. approximation to F Q with respect to G^. By 2.3 and 
Brouwer’s theorem, each f^ has a fixed point x . Hence 



x* = T X .y*’ 1 = I. X .f^ 5 ^ for X . > 0, £ . . T X = 1 (*) 

^i£N 0 k,r Z ' 1 ^ N 0 k » x k,i - kjl 



with <y k, °,. . . ,y k,n ) € G k , f* ,x € T Q (y^ iX ) for each i € N Q . 

As k -► °°, x , all the X .’s and f ’ ’s remain in compact sets C, [0,1] 
k ,i 

and C. Thus there is a subsequence (without loss of generality, the whole sequence) 
k k x d. 

along which x -»■ x 5 *, X^ ^ ■+ X^ and f 5 -*■ f . Since mesh G^ -*■ 0, we have that 
k i 

y * + x* for each i £ N . Because F Q is u.s.c., we deduce from 1.3 that 
f 1 £ F (x*) for each i £ N. Taking limits in (*), we have x* = N X^f 1 , 



■ k 5 /“ r 



x. > o, r . x. 



1. 



But F (x*) is convex; hence x* lies in F (x*) and the result is proved. | [ 



2 ♦ 4 Remarks ♦ It is crucial to note that x cannot be an approximate fixed 

point of f^ generated by the algorithms of Chapter IV. Indeed, consider the fol- 
2 

lowing example where C = S . Let F(z) = D for all z close to x. Then however 
fine the triangulation G, there is a p.l. approximation f to F with respect to 
G so that a simplex containing x is completely-labelled by the standard labelling 
of Chapter IV. (The diagram shows how such a p.l. approximation f can be 
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2 

v 




constructed.) However, clearly x is not a fixed point. 

It is also important to notice the ease with which the mapping F was extended 
to F , compared to the difficulties using retractions in II. 2. 

We develop algorithms to find the fixed points of p.l. approximations later. 

Let us establish how close such points are to their images under the original mapping. 

2.5 Theorem . 

(a) Let C C R n be triangulated by G, with mesh 2 G <_ 6 . Let F : C + C* 

satisfy the condition: for all x £ C, y £ B(x,6), we have F(y) CB(F(x),e). 

Let f be a p.l. approximation to F w.r.t. G and x* be a fixed point of f. 
Then x* £ B(F(x*),e). 

(b) Let F: C + C* be as in (a). Pick c £ int C. Let F^x) = F(x) if 

x € C, 0 if x £ C, F (x) = {c} if x £ int C, 0 if x € int C, and 

F q : R n ■+ C* be defined by F Q = conv(F 1 U F 2 ). Let R n be triangulated by G 
with mesh 2 G <_ 5 , f be a p.l. approximation to F Q w.r.t. G, and x* be a 
fixed point of f. If B(c,y) C C C B(c,v) with y > z , then 
x* g B(F(x*), e + (5+e)(v+e)/y) . 

(c) Let C C R n be triangulated by G with mesh 2 G <_ 6. Assume dim C = n. 

Let f q : C -* C be continuously differentiable with derivative f^ satisfying 

|| f q ( x ) - fg(z)|| _< M||x - z|| 2 for all x,z e C. Let f be the p.l. approximation 

1 2 

to { f Q } w.r.t. G and x* a fixed point of f. Then ||f Q (x*) - x*|| 2 < - M5 . 

Proof . 

.(a) We have x* = X.y 1 = 1.^ X.ffy 1 ) with X. > 0, X. = 1 
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and (y°,...,y n ) £ G. (We assume without loss of generality that dim C = n.) Now 
Hy 1 - x* || <5, so fCy 1 ) £ F(y 1 ) CB(F(x*),e). Since B(F(x*),e) is convex, it 

contains all convex combinations of the fCy 1 )^, in particular, x*. 

(b) By the argument of (a) we have x* = Ac + (l-A)d for some 0 <_ A _< 1, 
where d £ B(F(x*),e) and X - 0 unless x* is within 6 of R n ^ int C. If 
X - 0 we are done; hence assume X > 0. Because x* - d = A(c-d), we have 



x* £ B(F(x*), e + A j| c - d|| 2 ). 

Now F(x*) CC, so d £ B(C,e). Let d r £ C fl B(d,e). Since B(c,u) C C, 

B( Ac + (l-A)d' ,Ay) CC. Also, ||x* - (Ac - (1-A)d , )|| = (1-A) j|d - d T || ^(l-A)e. 
Because x* is within 6 of R n ^ int C, we must have 6 >_ Ay - (l-A)e >_ Ap - e. 

Hence A <_ (S+e)/y and ||c - dj| 2 <_ ||c - d* \\^ + ||d’ - d|| 2 <_ v + e . 

(c) For any x,z £ C, we have f^(z) - fg(x) = ^0 + = 

/J [^(x + A(z-x)) - fg(x)](z-x)dA + fg(x)(z-x) . Hence 



Il f 0 (z) - f 0 (x) - fg(x)(z-x)|| 2 = ||/g Cfg(x + X(z-x) ) - fg(x)](z-x)dX|| 2 

i Iq l|fo (x + X(z ' x)) - f o Cx) ll 2 H z ■ X ll 2 dx 

l/o MX H Z - X ll 2 dx = f M ll z - x ll 2 - 



Now 



x* = L^r ^-y 1 = F.^ T A.f (y 1 ) 
L i€N Q i y Ll €N 0 i 0 vjr 



= W W*** 



= + I 



ie N o 



Since each summand has norm at most 



with x i i °> Z ieN x i - (y°»- • • ,y n > € g 

+ ^igN 0 x i (f o (yl) ‘ f o (x ’ V) ' f o (x * )(yl ‘ xA)) 

+ h &0 x i f ;< xft)(yi - x *> 

x i (f o (yi) ■ f o (xft) * - **>) 

— A^M6 , the result follows. | | 



V.3 Extensions of Kakutani’s Theorem 

3 . 1 Corollary (Eaves [9]). Let C CR n be compact and convex with c£ int C. 
Let F: C R n be a u.s.c. mapping satisfying c £ F(x) whenever x £ 3C. Then 
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F has a fixed point in C. 

Proof . Since F is u.s.c., F(C) and hence C f = conv(C U F(C)) is compact. 
Extend F to F^: C* C T * using c £ int C in the same way F was extended to 
Fq in 2.1. Fq has a fixed point x* by 2.1. If x* £ C, Fq(x ! ' { ) = {c} C C, a 
contradiction. Thus x* £ C and x* £ F (x*) = F(x*). Q 



The next result is useful in the computation of economic equilibria. 

3.2 Corollary . Let T n = aff(S n ) = {x € R n+1 |v T x = 1} and let F: S n •+ T n " 
be u.s.c. Suppose there exist f 1 £ T n - T n , i £ N Q , such that 

(i) x ^ sS x + f 1 € F(x); and 

(ii) f j < 0 for j 4 i. 

Then F has a fixed point in S n . 



Proof . Let S n = {x g T n |x > -v}. Define u.s.c. mappings F^ , i = 1,2,3, as 

follows. F^(x) = F(x) if x ( S n , 0 if x £ S° 'v S R ; F 2 (x) =0 if 

x £ rel int S n , {x} + conv^^i ^ 1} if x £ S U ^ int S n and I = {i|x^ = min_.x_.}; 

and F^(x) = 0 if x £ rel int S n , {(n+1) ^v} if x £ 3S n . Finally, let 

F q = conv(F U F 2 U F 3 ) : S n -* T n5 \ Then F Q is u.s.c. and (n+1) 1 v € F Q (x) for 

o, n , . ^n 

x £ 9S . By 3.1, F Q has a fixed point x* m S . 

If x* £ S n , let I = { i | x* = min^x^*} . Let z be any member of F Q (x*) C 

conv{(n+l) ^v, {x* + f 1 | i £ I}}. We show that z j > x j an ^ hence 

x* 4 z. 



If z = (n+1) 1 v, z^ > 0 > x*. If z = x* if 1 , i € I> then 



L 



j£l j 



I j€ I (x* + f}) = l. & X* - l f* > I jeI x* by (ii); note that 

-1 U “1 
I 4 N q since x* 4 (n+1) v. Since the inequality holds for z = (n+1) v and 

z = x* + f 1 , i £ I, it holds for any z in their convex hull F Q (x*). 

We conclude that x* € S n so that x* £ F Q (x*) = F ^ x5 ‘^‘ CH 



For n = 2, a typical choice of f^’s is f° = (l,-j,-^-) T , f 1 = ( - y, 1 , - j) T , 
2 1 1 T 

f = (--,-— ,1) . Then the hypothesis of 3.2 requires that points on the boundary 
2 

of S are mapped as shown below. 
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3 . 3 Corollary . Let F: R n ■+ R n " be u.s.c. Suppose there is x° £ R n and 
y > 0 such that whenever x £ B(x^,y) there is some w £ R with w (x -x) > 0 

T 0 

and w (f-x) > 0 for all f £ F(x). Then F has a fixed point in B(x ,y). 

Proof. Let C be the compact convex set B(x^,2y). Let F^ be the restric- 
tion of F to C, F 2 be 0 on int C and {x°} on 3C, and F Q be 

n * 0 

conv(F^ UF ). Then F Q is a u.s.c. mapping from C to R with x £ F Q (x) 
whenever x £ 3C. By 3.1, F Q has a fixed point x*. If x* £ B(x°,y), F Q (x*) = 
F(x*) and the result is proved. 

Suppose x* £ B(x°,y). Then x* £ F Q (x*) C conv{x° ,F(x*) } . Thus x* = Xx° + 
(l-X)f for some X £ [0,1], f £ F(x*). But then there is w £ R n with 
w T (x°-x*) > 0 and w T (f-x*) > 0; hence w T (x*-x*) >0, a contradiction. Q 

Notice that the result also holds if the condition is required only for 
x £ 3B(x°,y). Since the condition is equivalent to x - x° £ U^ >0 X(F(x) - {x}), 
the result is then an extension of that of 1.5.2. 

For practical computation we need a slight strengthening of 3.3. The following 
result is due to Merrill [48,49]. 

3 . 4 Corollary ♦ Let F: R n •+ R n " be u.s.c. Suppose there are x° £ R n , y > 0 
and 6 > 0 such that whenever x £ B(x°,y), f£F(x) and z £ B(x,6), 

(f-x) T (x°-z) > 0. Then F has a fixed point in B(x°,y). 
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Proof . The condition implies that w = x 
3.3. □ 



z satisfies the hypotheses of 



Again we can strengthen 3.4 by requiring the condition to hold only when 
x £ 8B(x°,y ), but in this case the result is unsuitable for computation. If the 
condition need only hold for x £ B(x°,y+26) ^ B(x°,y), then a result is obtained 
which is appropriate for computation; on the other hand it is not clear whether the 
applicability of the result is appreciably extended for practical problems. 

We will develop algorithms based on 3.1 and 3.4. For those based on 3.1, we 
do not need C explicitly. We must have available a half space containing C and 
c £ int C, and, for any x £ R n , we must be able to determine whether or not x 
lies in C. The algorithms based on 3.4 require only a knowledge of 6. Fortu- 
nately, for many problems the condition of 3.5 holds for any 6 > 0 with appropriate 
x° and y . Clearly, if F(R n ) is compact (as it is, for example, if F is the 
natural extension to R n of a mapping satisfying the conditions of 3.1), then this 
stronger condition holds. For all algorithms, as indicated by the proof of 2.1, 
we need only find a single member of F(x) for each x £ R n . 




CHAPTER VI: APPLICATIONS OF KAKUTANI’S THEOREM AND ITS EXTENSIONS 



We now show how the ability to find approximate fixed points of upper semi- 
continuous mappings enables one to solve a wide range of important problems, without 
the limitations of Chapter II. Again we point out that constructing an appropriate 
mapping is relatively straightforward. Indeed, for any given problem it is easier 
than constructing a continuous function, since u.s.c. mappings can be patched 
together so easily. 

VI. 1 Equilibria in Finite Economies with Production . Consider the model of an 
exchange economy as in II.l. The demand side of the economy we will consider is the 
same except for the following important relaxation. For any p £ S n , we assume that 
the ith consumer has a corresponding set D 1 (p) of possible demands. D 1 is assumed 

to be a u.s.c. mapping from S n to R^ +1 (see exercise 5.1); set D(p) = 

r*m i n n+l^' 

2, i=1 D (p). The aggregate demand correspondence D: S -*• R + is then u.s.c.; 

assume that it satisfies Walras' Law: for all p £ S n , d £ D(p), p T d = p^w. 

On the production side we restrict ourselves to a simple model. Let B = [-I,C] 

be an (n+l)x(k+l) matrix. For 0 <_ j _< k, the jth column B_. of B represents 

an activity; if activity j is used at level 1, | b „ | units of commodity i are 

supplied as output (if b^_. >_ 0) or required as input (if b^ < 0). Note that the 

first n+1 columns of B correspond to an assumption of free disposal. 

The set of all production vectors is {Bz|z >_ 0). The model can be generalized 

to allow more general production sets, but they must at least be convex. We make 

the reasonable assumption that {Bz|z _> 0} f) R^ +1 = { 0 } ; production is bounded. 

Equivalently, one assumes that there is q > 0 with qB < 0. 

We can now define an equilibrium as a price vector and a vector of activity 

levels, such that the sum of the initial resources and the production vector is a 

member of the demand set for those prices. (Since free disposal is included in the 

production possibilities, we need not worry about excess supply.) By permitting 

only activities that maximize profit we ensure that producers have an incentive to 



use them. 
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1.1 Definition . An equilibrium is a pair (p*,z*) £ S n x R^ +1 such that 
w + Bz* £ D(p*) and p ssT B <_ 0. 

Using Walras' Law we find that p* w + p" Bz* = p* w, and hence p rt Bz* = 0, 
p*^B < 0. The interpretation is that no activity makes a positive profit; those 
activities used (at positive level) must make zero profit. In particular, if a free 
disposal activity is used, the price of the corresponding commodity is zero. 

We shall prove that an equilibrium exists (if w > 0) by constructing a u.s.c. 
mapping satisfying the hypotheses of V.3.2. Scarf [58] does not require w > 0 and 
does not need the vector q above, as we do. However, our mapping allows the use 
of more sophisticated algorithms for fixed-point problems. 

1.2 Theorem. If w > 0 then an equilibrium exists. 

Proof. Let p £ S° be given. If p is not an equilibrium price vector, we 
wish to modify it suitably. First, if any p j _ = 0 , we may want to increase p. to 
avoid leaving S n . Secondly, if p > 0 but tt* = max^p^ >0, we may want to 
adjust p to decrease the profit of the most profitable activities. Finally, if 
p > 0 and tt* < 0, we adjust p to decrease demand. 

Define the following mappings: 



E^P) 



conv{v^|p^ =0} if p € 3S ; 
0 otherwise . 



E 2 (p) 



T T 

conv{-B Ip B = tt* = max p B 
r r s s 



if tt* > 0; 



otherwise . 



e 3 ( p) 



D(p) if it* 0; 



0 otherwise . 



Clearly, 

mappings such as 



is u.s.c. for 




i = 1,2,3; for example, 
if p T B r > 0, p T (B r -B s )>_0 

otherwise 



is the convex hull of 
which are u.s.c. by 
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V. 1.4(e). Now let E = conv^ U E 2 U E^. Then E: S n + R n+1 ' k is u.s.c. We would 
like to send p to {p} + E(p), but unfortunately this set may not even lie in 
T n = {x £ R n+ 1 |v T x =1}. We must modify E(p) so that e £ E’(p) => e £ T° - T n . 

For any e £ E(p), proceed as follows. 

Let e’(e) have coordinates e^, i £ N Q , defined by 

e! = q.e. (q T w/q T e) - q.w. . 

1 l l ^ l 

To show that e’ is well-defined and e’(*) is a continuous function, we need 

T i T 

only show that q e is bounded away from zero. But for e=v, q e = q. > 0. 

For e = -B , q^e = -q T B > min (-q T B ) > 0. For e = d f D(p), q^e > (min.q. )v T e 
r ^ r — s s — i n i 

T T 

> (min.q.)p e = (min.q.)p w > (min.q. )(min.w. ) > 0. Since E(p) = 

— ii ii — ii 33 

T T 

conv(E (p) U E (p) U E (p)), we have q e >_ min {min.q. , min (-q B ), (min.q. )x 

-L Z 11 S S 11 

(min.w.)} > 0 . 

3 D 

Thus we may set E'(p) = {e’(e)|e £ E (p ) } . By V. 1.4(b), E’ is u.s.c. Note 

that e'(e) lies in T n - T D for any e with q^e i 0 ; it is an easy exercise 

(5.2) to show that E’(p) is convex. Thus E' : S n -► (T n - T n )*, and F: S n + T^' 
is u.s.c. with F(p ) = {p} + E'(p). 

Now if p £ S?, v 1 £ E(p), so e^v 1 ) £ E (p). For each i £ N^ define 

f 1 = e ' ( v 1 ) = (-q Q w 0 , . . . ,-q^w^ + q T w, . . . i~<I n w n ) T - Thus for p £ S?, p + f 1 £ F(p). 
Thus the hypotheses of V.3.2 are satisfied, and F has a fixed point p* £ S n . 

Then 0 € E'(p*), that is, there is e £ E(p*) with e = Xw for some positive 
scalar X. From the definition of E we have 

Xw = £. y . v^ + v (-B ) + pd (*) 

ii^N Q i ^r =0 r r 

with y. , v , p > 0, y y. + 7 v + p = 1. In addition, y. = 0 unless p? = 0, 

l r — ^ i ^ r l r i 

.'.T a T 

v =0 unless p 5 ' B = max p 54 B = tt* > 0, p = 0 unless it 5 ’* < 0, and d G D(p*). 

r r r s r s — - 

A t 

Premultiplying (*) by p“ yields 

Xp* T w = 0 + (l - Ott* + pp* T d. (**) 
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If tt* > 0 , p = 0. Then the left-hand side of (**) is positive and the right- 
hand side nonpositive, a contradiction. 

If ir* < 0, all = 0, and (**) becomes Ap**w = pp“ T d. Using Walras' Law, 
we obtain A = p. Dividing (*) by X gives 

w + h € N 0 = d- 

T 

Hence (p* , (p Q /A , . . . /A ,0 , . . . ,0) ) is an equilibrium; the only activities 

possibly used are those of free disposal. 

Finally, if it* = 0, then again Ap'^w = pp*^d and X - p. Since v 1 = -B^ , 

i € N q , and pV = 0 implies that activity makes the maximum profit zero, we 

can incorporate all the p.'s that are non- zero in the v 's. Then from (*) we 
i r 

have 



w + J* (v /X)B = d. 

^r=0 r r 

Hence (p*, ( v^/A , . . . ,v^/X) ) is an equilibrium and the theorem is proved. ] [ 

VI. 2 Unconstrained Optimization . The next three sections follow Merrill's analysis 

[48,49] closely. We restrict ourselves to the simplest case. 

Let 0 : R n -*■ R be convex. We write lev 9 for {x f R n j9(x) < a} and D9 

a 1 — 

for the subdifferential mapping of V.l. Then with 0 finite everywhere, we con- 

n n * 

elude that 9 is continuous, lev 9 is closed and convex, and D0 : R -*■ R is 

a 

u.s.c. (Rockafellar [51], theorems 10.1.1, 4.6, 23.4 and 24.5.1). Also, if l ev a e 
is nonempty and bounded for some a, it is bounded for all a (corollary 8.7.1 of 
Rockafellar [51]). 

Define : R n -> R n by F^(x) = (x) - D0(x). Let 9(c) < a, and define 

F'(x) to be F_ (x) on lev 0, 0 otherwise, F"(x) to be {c} on R n ^ int lev 9, 
la a 

0 otherwise, and F ^ = conv(F' U F"). 

2 . 1 Lemma . The set of fixed points of F^ is the set of minimizers of 0 , 



for i = 1,2. 
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Proof . The result is clear for For F 2 we need only verify that F 2 has 

no fixed points on 9 lev a 9 * But x € 9 lev a 9 ’ x ^ ^(x), then x = + 

(l-A)(x-d), for A £ [0,1] and d £ D6(x). Then d = A(l-A) *(c-x) and 0(c) > 
0(x), a contradiction. □ 

2.2 Theorem . Suppose that 0 has a bounded nonempty level set. Then F^ 
restricted to C = lev^O satisfies the hypothesis of V.3.1 and F^ satisfies the 
hypothesis of V.3.4. Indeed, for any x 9 £ R n and 6 > 0 there is y > 0 satis- 
fying the condition of V.3.4. 



Proof . Since C is compact, the first part is clear. Now let 
3 = sup{0(x)|x £ B(x 9 ,5)} < Then lev^0 is bounded; choose y so that lev^0 S 
B(x°,y). Let x £ B(x°,y), f £ F^Cx) and z £ B(x,6). Then x - f £ D0(x) and 
hence (x-f) T (x°-z) = (x-f ) T ( (x°-z+x) - x) <_ 0(x°-z+x) - 0(x). But x° - z + x £ 
B(x°,<5); thus (x-f) T (x°-z) <3-3=0. Q 



If a fixed point of a p.l. approximation to F^ is found, a lower bound on 0 

is naturally generated. If x* is such a point, x* = 7 . N A y 1 = A (y 1 -d 1 ), 

lfcW 0 1 0 1 

where A. > 0, T . A. = 1 and d 1 £ DOCy 1 ) for each i £ N_. For every x £ R 
i — 1 ^0 1 ° 

and each i £ N^, 0(x) 0 ( y 1 ) + d 1 ^(x-y 1 ); hence 



V (yi) + W x i diT(x - yi) 



= x i e(yl) + W x i dl (x *- yl) (A) 



since £ A.d^ = 0. If diam(y^ , . . . ,y n ) is small, the last term in (*) is likely 
to be small; the first term is at least 0(x*), showing x* to be close to 
minimizing 0 . 

Merrill [48,49] discusses extensions of this application, including the case 
where 0 is not convex. 



VI. 3 Constrained Optimization. Consider the problem: minimize 0(x) subject to 
g k (x) < 0, k = 1,2 , . . . ,m, where 0 and each g k is convex and takes R n to R. 
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By defining ^ ( x ) to be max^g^Cx), we reduce this problem to 

min{9(x) | ip ( x ) <_ 0} , (P) 

where 0 and ^ are convex. If each g^ is differentiable, then Dtp(x) = 
conv{Vg^(x) |g^(x) = iKx)}. 

A natural modification of any x £ R n that does not solve (P) motivates the 
following definitions. Let 

{ {x} - D9(x) if x £ lev^, 

0 otherwise ; 



H 2 (x) 



0 if x £ int lev^, 

x} - D ip ( x ) otherwise; 



and ? 1 = conv(H 1 U H 2 ): R n + R n '\ 

If there is a c £ R n with 4>(c) < 0, we define F^ as follows. Let 9(c) 

< a, i|/(c) < 0 < 3, and C = lev 9 fl 1 ev a ip. Define F T (x) to be F (x) if x € C 

a p l 

and F"(x) to be 0 if x £ int C, {c} if x € 3C. Finally, let F 2 = 
conv(F' U F M ): C + R n '\ 

We write inf ip to denote inf{ij;(x)|x £ R n }. The relationship between F^, 

F 2 , and (P) is shown by 

3 . 1 Lemma . Any fixed point of F^ lies in lev Q ijj if this set is nonempty. 

If inf ip < 0, then the set of optimal solutions to (P) is the set of fixed points 
of F. , i = 1 or 2 . 

l 

Proof . Let x* £ F^Cx*) ^ lev Q ^. Then 0 £ D^(x*) and ip takes on its 
minimum at x*. Thus inf ip = \p(x' e ) > 0 and lev^ is empty. 

Now assume inf ip < 0. Let x* be an optimal solution to (P). It follows 
(see Rockafellar [51], corollary 28.2.1) that x* minimizes 9 + Xip -for some 
X >_ 0 with X = 0 unless i|>(x*) = 0. Thus 0 = (1+X) ^d + X(l+X) ^e with 
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d £ D0(x*), e £ DiJ>(x*), and x* £ F^tx 5 *). Since 0(x*) <_ 6(c) < a and i p( x* ) <_ 

0 < 3, x* £ int C and x* £ F^Cx*) also. 

Conversely, let x* be a fixed point of F^. If ip(x*) < 0, 0 £ D9(x*) and 
hence x* solves (P). If \p(x t,{ ) = 0, 0 = Xd + (l-X)e for some X £ [0,1], 
d £ D0(x*) and e £ Dip(x *). If X = 0, x* minimizes a contradiction. 

Dividing by X > 0 we obtain d + X "''(l-X)e = 0, and x* minimizes 9 + X 1 — X )t(i . 
Thus x* solves (P). 

Let x* be a fixed point of F^. Clearly, x* £ C, and if x* £ int C, x* 

is a fixed point of F^ and hence an optimal solution to (P). Let x* £ 3C. If 

x* £ lev ij/, we obtain a contradiction as in 2.1. Thus <K X *) <_ 0 < $ and 
0(x*) = a. If < 0, then 0 £ conv(D0(x*) U {x-c}) and a contradiction 

arises as in 2.1. If ij>(x*) = 0, then 0 £ conv(Dx(x*) U {x-c}) with x = ^9 + 
(l-X)i^, some X £ [0,1]; again a contradiction results from 0(x*) = a > 0(c), 
i|>(x*) = 0 > i|;(c) and hence x(**) > x(°)- EU 

If lev 0 D lev J is nonempty and bounded for some a, 3, it is bounded for 
a p 

all a, 3. (Probably the easiest way to see this is by considering functions of the 
form cf>(x) = max{X(0(x) - oc), y(ip(x) - 3)}, and applying corollary 8.7.1 of 
Rockafellar [51].) This condition is sufficient for the existence of fixed points 
of F^ and F^. More precisely, we have 

3.2 Theorem . If ip has a nonempty bounded level set, then for any x° £ R n , 

6 > 0 there is a y > 0 satisfying V.3.4 for F^. If inf ip < 0 and lev^0 fi 
lev^ is nonempty and bounded for some a, 3, then satisfies the conditions 

of V.3.1, and there are x° £ R n , 6 > 0, and y > 0 satisfying V.3.4 for F^. 

Proof . If ijj has a nonempty bounded level set, then lev Q i^ is bounded and 

F^(x) = H^(x) for x sufficiently far from 0. The first conclusion then follows 

as in 2.2. The result for F^ is clear. If inf ^ < 0, let ij/(x°) < 0 and 6 

be such that 3 = sup{if»(x)|x £ B(x°,6)} < 0 and let a = sup{0(x)|x £ B(x°,5)}. 

Choose y so that lev 0 D lev„i|; CB(x°,y). The conclusion follows as in 2.2. I I 
a 3 — 1 — 1 



For computation, the first part of 3.2 is much more satisfactory than the last, 
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for we choose 5 arbitrarily. If it is known that there is a feasible solution to 
0(x) < a, ip(x) < 0, then (P) is equivalent to: minimize 0(x) subject to 
iKx) 0, 0(x) <_ a, or minimize 0(x) subject to \p'(x) £ 0 (P’), where 
i|>'(x) = max{i|>(x), 0(x) - a}. If lev^ T 0 fi lev^ij; is bounded and nonempty for 
some a 1 , 3', then t/i ' has a bounded nonempty level set. Thus basing on (P f ) 

rather than (P) allows one to choose 6 arbitrarily. 



VI. 4 The Nonlinear Complementarity Problem 

Recall that the NLCP associated with a continuous function g: -► R n is to 

find x* _> 0 with g(x*) > 0 and x”^g(x*) =0. We describe how Merrill [48,49] 
converts the NLCP into a fixed-point problem and give sufficient conditions for the 
resultant mapping to satisfy the hypothesis of V.3.4. 



Let ip : R -* R be the convex function with i|j(x) = max. __(-x.); note that 

l£N 1 

lev^ = R^. Define F^(x) to be 0 if x £ int R^ and {x} - Dt|j(x) otherwise, 

and F 2 (x) to be (x - g(x)} if x £ R^ and 0 otherwise. Let F = convCF^ U F 2 ): 
n* 

R -> R . 



4 . 1 Lemma ♦ x is a fixed point of F iff x solves the NLCP associated with 

g- 

Proof . Let x be a fixed point of F. If x £ R^, 0 £ Df(x) and x mini- 
mizes ip, a contradiction. If x £ int R^, g(x) = 0 and x solves the NLCP. 

If x £ 9R^, set I = (i £N|x^ = 0}. Then Dip ( x ) = conv{-u 1 |i £ 1} and hence 
g^(x) = 0 for i £ I and g^(x) >_ 0 for i £ I. Thus x solves the NLCP, 

Conversely, given a solution x to the NLCP, set l={i|x^=0}. If 1=0, 
g(x) = 0 and x £ F 2 (x) = F(x). Otherwise, -g(x) £ conv{-u 1 |i £ 1} = D^(x) and 
x = i<x + g ( x ) ) + j(x - g ( x ) ) £ conv(F^(x) U F 2 (x)) = F(x). Q 

A sufficient condition for F to satisfy the hypothesis of V.3.4 (and thus 
for the NLCP to have a solution) is given by 

4.2 Theorem. Suppose there exist a > 0 and v > 0 such that whenever 

x (= R^ ^ B(0,v), x T g(x) > a ||x|| 2 ||g(x) || 2 - Then for any 6 > 0 there exist x° £ R n 
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and y > 0 satisfying the hypothesis of V.3.4 for F. 

Proof . Choose x° > 6u and y > max{v + ||x°|| 2 , (6 + ||x^|| 2 )/a}. Note that 
if;(x) < 0 for any x £ B(x°,6). Now let x £ B(x°,y), f £ F(x), and z £ B(x,6). 
If ip( x) > 0, then -(f-x) £ Dip(x) . Thus, as in 2.2, we find (f-x) T (x°-z) _> 

^(x) - ^(x^+x-z) > 0. 

If iKx) <0, x £ R^ and x £ B(0 , v) . In this case f = x - g(x) and 

(f-x) T (x°-z) = -g(x) T (x°-z) = g(x) T x + g(x) T (-x°-x+z) > ot ||x |) 2 ||g(x) || 2 - 

||g(x)|| 2 ||-x°-x+z|| ■> ||g(x) || 2 (a||x|| 2 - i|x°|| 2 - 6) > 0, since g(x) cannot be 0. 

Finally, if i|>(x) = 0, a combination of the cases above yields 
(f-x) T (x°-z) >0. | | 

VI. 5 Exercises . 

5 . 1. Under the assumptions of II. 4.1, show that D is a u.s.c. mapping on 

t" = {p e s n | P > o}. 

5.2. Show that E’(p) in the proof of 1.2 is convex. 

5.3. Let 0 ,ip : R n -► R be convex with lev^B (1 lev^ f i|; nonempty and bounded for 

some a', 3'. Show that lev 9 fl lev J is bounded for all a, 3. 

a p 

5.4 . Let 0 : R n R be convex with a bounded nonempty level set. Pick x\x 2 £ R n 

with OCx 1 ) ^ 0(x 2 ). Using only x 1 , ©(x 1 ), and some d 1 £ DBCx 1 ) for i = 1,2, 

show how to obtain: a so that C = lev^O is compact; some c £ int C; and a 

halfspace containing C. 

5.5 . Let 6,ip: R n -+ R be convex with lev a i 9 fi lev gt^ nonempty and bounded for 

some a', 3'* Pick x\x 2 £ R n with ^(x 1 ) < 0 and 0(x'*') i 0(x 2 ) if 

\ff(x 2 ) <_0. Using only x 1 , ©(x 1 ), ^(x 1 ), some d 1 £ DOCx 1 ), and some 

e £ Di^Cx 1 ), i = 1,2, show how to obtain: 3 > 0 and a so that 

C = lev 0 fl lev.ip is compact; some c £ int C; and a halfspace containing C. 
a 3 

5 .6. Suppose the condition of theorem II. 3.1 holds for some w > 0. Use V.3.1 to 

prove that the NLCP associated with g has a solution in C . (Hint: w £ int C^.) 




CHAPTER VII: EAVES’ FIRST ALGORITHM 



The first algorithms for u.s.c. mappings were developed by Hansen and Scarf 
[29]. However, it is more natural for us to use triangulations rather than primi- 
tive sets, and hence we will describe Eaves' algorithm [10,12]. Section 1 intro- 
duces the crucial notions of vector labelling and complete simplices. Section 2 
describes how the problem to be solved, based on V.3.1, is set up. In Section 3 we 
define a graph T similar to the graph of Chapter III. Eaves’ algorithm is 

presented in Section 4. 



VII. 1 Vector Labelling and Complete Simplices . Let A be an affine subspace in 
R m , and suppose dim A = n. Let F: A -► A* be u.s.c., G be a triangulation of 
A, and f be a piecewise-linear approximation to F with respect to G. Let S 
be the linear subspace A-A of dimension n. Define l: A -*■ S by £(x) = f(x)-x. 
We call l the vector labelling for f (or a vector labelling for F with respect 
to G). Often, l is taken to be a function on the vertices of G; however, 

extending l to R n shows that fixed points of f are zeroes of i . 

0 k k 

If a = (y , . . . ,y ) £ G , then the (m+l)x(k+l) matrix 

r i ... i i 



L = 



.«Cy°)...«.(y k ). 



is called the label matrix of a. 



Let t° be the unit vector (1,0,...,0) T in R m+1 



fixed points of i by means of 



Then we can identify 



1 . 1 Lemma . Let A, F, G, f, and J l be as above, and let a = (y°,...,y n ) £ 
G. Then there is a fixed point of f in a if and only if there is a solution to 



L w = t°, w £ R n+1 . 
a + 



(*) 



If w solves (*), then x* = w.y 1 is a fixed point of f 






Proof. For any x = 7. .. w.y 1 a a we have 
J i 



[ i@I v (yl) = Hx) = f(x) ' x 



(t> 
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If x* is a fixed point of f in cr , 

„ "-y 1 ; an< ^ by (t) w is a solution to 

1 ^ N 0 1 

Now let w solve (*) and set x* = J. 

lfcN 

from the first equation in (*); thus x* £ a. 



then we can find w £ S n with x* = 
(*). 

i T 

w.y . Then w > 0 and v w = 1, 

0 1 

Now from (t) we find f(x*) = x*. Q 



Motivated by 1.1, we state 



1.2 Definition . A simplex a of G is complete if there is a solution to 

(*). 



Vector labelling can be regarded as an extension of integer labelling, and 
complete simplices of completely-labelled simplices, as follows. Let g: S n ■+ S n 
be continuous. Suppose we use the integer label min{i|g^(y) - = min^(g^(y )-y^) } 

for y. This labelling is admissible if g has no fixed points on 3S n . Now let 
A = aff(S n ). For x £ rel int S n , let I(x) = (i|g^,(x) - x^ = min^g^Cx)-^) } . 

For x £ 3S n , let I(x) = {i|g^(x) - x^ = min^(g^(x)-x^) } U {i|x^ = max^x^); and 
for x £ A ^ S n , let I(x) = { i J x^ = max^x^}. 

For each i £ N^ let r 1 be the ith column of the (n+l)x(n+l) matrix 
~-l 0 • • - 0, +1" 

+1 ‘ € * . ‘ 0 

0 ' , * 



“ . Define F: A -► A* by F(x) = {x} + convlr^i £ I(x)}. 



‘ # . -1 0 

o • . . 0 +1 -1 _ 

For each i £ N , {x £ A | i £ I(x)} is closed; hence F is u.s.c. It is easy 
to check that the fixed points of F are exactly those of g. 

Now let G be a triangulation of A. For y £ G° choose f(y) = y + r 1 with 
i the least index in I(y), and let f be the corresponding p.l. approximation to 
F and l its vector labelling. Then if a £ G lies in S n , it is completely- 
labelled (with the integer labelling) iff it is complete (with the vector labelling 



l). 



The preceding discussion is based on Eaves [14]. 



VII. 2 Preliminaries . For simplicity we take the set A of Section 1 to have full 
dimension, i.e., m = n. The necessary changes for m > n are straightforward. 
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Suppose we have a compact convex set C C R n and c £ int C. A u.s.c. mapping 
Fq : C R n is given, with c £ Fq(x) whenever x £ 9C. We must have available c 
and some halfspace containing C; for simplicity assume C C {x £ R n |x^ > 0}. 

We now extend F Q to R n by setting F^(x) to be F Q (x) if x € C, 0 

otherwise; F^Cx) to be 0 if x £ C, { c } otherwise; and F = conv(F^ U F^): 
n n* 

R -* R . Then F is u.s.c. and has the same fixed points as F^ . We will find 
fixed points of piecewise linear approximations f to F with respect to G. Let 
G = 5K X or 5J . 

Pick a starting point x^ with x^ = 0. To enhance convergence, x^ should 
be close to a fixed point of F — but note that x^ £ C. In the next chapter we shall 
see how our starting point can be chosen wherever we wish. If necessary, perturb x^ 
so that it lies in some (n-1) -simplex of G n \ Define d = x^ - c. 

Let f be a p.l. approximation to F with respect to G, and l its vector 
labelling. A fixed point of f can be identified if we can find a complete simplex 
a £ G. 

To find c.l. simplices in Chapter III, we had to relax our requirements and 

search a. c.l. simplices. A similar relaxation is necessary here. Note that L^ 

is square ; hence the solution to L^w = v^ will possibly be unique — we do not have 

the required degree of freedom to produce paths . 

To see how to define almost-complete simplices, consider the analogy with 

integer labelling as in the last section. For the vector labelling there , a simplex 

a is a. c.l. iff there is a solution to Lw+(^)X=t^, w £ R n+ \ X > 0. Thus 

an + — 

r 

it is natural to define almost-complete simplices by including an artificial column. 

2.1 Definition . Let a = (y^ , . . . ,y^ ) £ G n ^ U G n . We call a almost- 
complete if there is a solution to 

Lw+(j)X=v°, w > 0 , X > 0 . ( * ) 

ad — — 

When a £ G n the linear system above has n+1 rows and n+2 columns . If 
the matrix [L^,(^)] has rank n+1, barring degeneracy, there will be two solu- 
tions to (*) with n+1 variables positive. We would like this tc be true in 
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general, but we cannot assume that degeneracy does not occur. In most linear 
programming problems this assumption is fairly harmless; if degeneracy occurs, the 
objective function usually acts to avoid cycling. In fixed-point algorithms, how- 
ever, the columns of are likely to be close if the grid size 6 is small; with 

roundoff error, degeneracy can easily occur. Since convergence depends critically 
on unique pivoting, we must devise a method of circumventing degeneracy. 

We do so by using lexicographic systems. A row vector e is lexicographically 
positive if e i 0 and e_. > 0, where j = min{i|e^ i 0}; we write e > 0. We 
say e is lexicographically nonnegative and write e>° if e >0 or e = 0. A 
matrix is lexicographically positive (nonnegative) if each of its rows is. 

Let A and B be real matrices of dimensions mx(m+l) and mxp, and suppose 
A and B have rank m. Then any solution X to the linear system AX = B, 

X > 0, has rank m; hence X has at most one zero row. If X has one zero row, 
the columns of A corresponding to the remaining rows of X are linearly inde- 
pendent and form a basis for R m . In this case we call X a basic feasible solu- 
tion to AX = B, X > 0. 

Noting that the null space of A has dimension 1 one can easily prove the 
following standard result of linear programming. 

2 . 2 Theorem . Let A and B be as above. If AX = B, X > 0 has a solution, 
it has exactly two basic feasible solutions. Q 

See Dantzig [8] for a further discussion of lexicographic rules. Note that one 
basic feasible solution can be obtained from the other by a linear programming pivot 
step . 

A natural way to obtain a right-hand side of rank n+1 in our systems is to 
replace v° by the identity matrix I of order n+1. 

2.3 Definition . Let a = <y° , . . . ,y k ) € G n_1 U G n . We call a very complete 
(v.c.) if there is a solution to 

L W = I, W > 0, 

a — 



(i) 
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and very almost complete (v.a.c.) if there is a solution to 

L W + (be = 1, W > 0, e > 0. (ii) 

ad — — 

VII. 3 The Graph r . 

3.1 Definition . The nodes of T are v.a.c. n-simplices of G lying in 

(x € R n |* > 0} and v.a.c. (n-l)-simplices of G + lying in {x £ R n |x^ = 0}. Two 

nodes of T are adjacent if one is a face of the other or if they share a v.a.c. 

facet. 

Let us first examine the degrees of nodes of V. If t is a v.a.c. (n-1)- 

simplex and a node of T, t is a facet of just one simplex of G in 

{x £ R n |x^ 0}; and clearly this simplex is v.a.c. Since no facet of t can be 
v.a.c., t has degree one. If a is an n-simplex of G and a node of T, we 

distinguish the cases where a is v.c. and where it is not. 

If a is v.c. , the system (ii) of 2.3 has two basic feasible solutions, on6 
of which is the solution of 2.3(i). Thus a has just one v.a.c. facet t. Whether 
t lies in {x ^ R n |x 1 = 0} (and hence is a node of T) or not (in which case, x 

is a facet of another node a' of T), we find that a has degree 1. 

If a is not v.c., then the two basic feasible solutions to 2.3(ii) show that 
a has exactly two v.a.c. facets. Hence in this case a has degree 2. 

We now see that T shares many desirable features with of Chapter IV. To 

obtain an algorithm, we would like to know that T has only a finite number of nodes 
and has a natural starting node. 

Let Tq be the (n-1) -simplex of G n 1 containing x°. 

3 . 2 Lemma . 

(a) T has only a finite number of nodes. 

(b) is the only (n-l)-simplex that is a node of T. 

Proof . 

(a). Let a be a node of T. If c has a vertex in C, then a CB(C,<5/nj, 



which is compact. Thus there are only a finite number of such simplices. Suppose 
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that no vertex of cr lies in C. Then L 

a 



where a 



<y° 



k \ 

»y /• 



L. c-y . . .c-y 

Since a is v.a.c. , it is almost complete and there exist w , . . . ,w^ 0 and 

X > 0 with y k w. + X = 1 and J k . w.Cc-y 1 ) + Xd = 0. We know that p = 

— L i=0 l L i=0 l J 

£ . w. > 0, for otherwise X = 1 and d = 0, a contradiction. Let y. = w./p 

^i=0 i * li 

and v = X/p > 0. Then 



rk i , 

I i= o = c + vd, 

i.e., c + vd £ a. Because ff C {x ^ R n |x^ > 0}, v <_ 1; thus a meets C* = 

{c + vd|o < v < 1} and a C B(C* ,6/n) . Since the latter set is compact, there is 

only a finite number of nodes of T with no vertices in C. The conclusion follows. 

(b) Any (n-l)-simplex t that is a node of T lies in {x € R n |x^ = 0} and 

has no vertex in C. As in (a), x must contain a point c + vd, 0 <_ v <_ 1; and 

since its first coordinate must be zero, t must contain c + d = x*“\ Hence 

T = V 

0 v* 1 

We must now show that is v.a.c. Let x = 2 w^y > where = 

i _ mi m -i m 

(y , . . . , y ). Then w > 0 and u w = — . We show that (w ,-i-) solves (*) of 2.1. 

1 T 1 

Indeed, w > 0 and — > 0 with u w + — = 1, while 



V (yl) + i d = L € n H i (c ' yl) + 1 d = I 



C - y X° + i d = °. 



Hence is almost complete. Also, w and -j suffice as the first columns of W 

and e in 2.3(ii) and are strictly positive. If L = [L q(^)] is nonsingular, we 
can obtain a solution to 2 . 3 ( ii ) . But certainly y\...,y n , c-d are affinely 
independent; hence so are d, c-y\ . . . ,c-y n , proving L is nonsingular. | | 

Combining our analysis of the degrees of nodes of T with 3.2, we obtain 



3.3 Theorem. Each connected component of the finite graph T is either a 
simple circuit or a simple path each of whose endpoints is either or a v.c. 
simplex. □ 
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VI 1. 4 The Algorithm . 

4.1 . We trace the simple path in Y with as one endpoint. 

Step 0 . Let a be the unique n-simplex of G in {x ^ R n |x 1 _> 0} with t q 

as a facet. Let y + be the vertex of that is not a vertex of t q . Form the 

linear system (as in 3.2) showing v.a.c. Set m = 1. 

Step 1 . Calculate the vector label &(y + ) of y + . Introduce the column 

( \ ) into the basis of the current lexicographic linear system. If the arti- 

Hy ) 

ficial column (^) drops from the basis, STOP — a is v.c. Otherwise, a unique 

d n 

column ( _ ) drops, where y is a vertex of a . 

*(y ) m 

Step 2 . Find the n-simplex a m+1 ° f G containing all vertices of 
except y . Let y + be the new vertex. Set m m + 1 and return to Step 1. 

4.2 Example (n = 2). Let C = {x £ R 2 x i ,x 2 - with c = € 

int C. Let ^ = (x £ C|x 1 ,x 2 <_ 2^} , C 2 = (x £ C|x x < 2y< x^ and C 3 = 

{x € C|x x > 2i}. Let d 1 = (^) T , d 2 = and d 3 = (-^,0) T . For 

i = 1,2,3; set F^(x) = {x + d 1 } if x € and 0 otherwise; and set F 4 (x) = 

4 

{c} if x £ 9C, 0 otherwise. Finally, let F = conv(lL =1 Fj. One can easily 

1 1 T 

check that (2— ,2—) is the unique fixed point of F. 

Pick x° = (0,l|-) T so that d = x° - c = (-1,0) T . We will triangulate R 2 

0 2 
with J 1< Then x = ((0,1), (0,2)). The relevant portion of R , together with 

the labels of the vertices shown by arrows, is indicated below. 




The initial simplex is a = <(1,1) T , (0,1) T , (0,2) T ) = (y^y^y ) = j 1 (y° ,tt ,s ) 
with tt = (1,2) and s = (-1,1). We have t° = (y\y 2 ). Let B = [L Q , (^)] = 



r 




1 1 


111 




=■ - 1 






4 4 




-1 


1 1 


1 — 1 
i — 1 
' — \ 


be the current basis. Then B = 


- - -1 

4 4 


1 I o 




1 .1 o 


2 2 ° 




- 2 2 



1 2 

first and second rows correspond to y and y and the third to the artificial 



column (^) (hence the notation M a"). Since the rows of B are lexicographi- 
cally positive, is v.a.c. The iterations now proceed as follows. 





?*) HJ I N ?“) >0 I 



H | CO <H |cO 04 |C0 



1^- h|j- H |cn 



■ 



i h|co h|j h o ol n|in to|m h|io h o o o rH o H |co H |co oj |co 



— I | 04 r- 1 | J H | 04 00 | J- 
H 04 H 04 



04 H 00 Ol 04 04 CO CO CO 



004 0) 004 CO 04040) Ol Ol 0) 0401 CO 



04 p) I CN CO f\| Ol CO Ol 



the artificial column drops from the basis, a is very complete. From the first column of the basis 



CHAPTER VIII: MERRILL'S ALGORITHM 



VIII.l An Extra Dimension . We will call the algorithms we have developed so far 
first-generation algorithms. They all suffer from computational inefficiency, caused 
by the following two features : 

(i) we must choose a particular grid size with which to work throughout the 
algorithm; and 

(ii) we must start outside the region of interest. 

If the grid size is large, we obtain a poor approximation, while if it is small, 
we spend an inordinate amount of time moving through simplices .r from the solution. 

One possible way to circumvent these difficulties is to refine the triangula- 
tion after finding a c.l. or v.c. simplex. The idea is illustrated with integer 
labels below: 

11 11 





Unfortunately, the refined triangulation within has simplices of as large 

diameter as itself, and this characteristic cannot be eliminated without 

allowing escape from a^. 

We will discuss two ways to avoid the problems of first-generation algorithms. 
The first method is that of Merrill and circumvents feature (ii) above; the algo- 
rithm can start anywhere. Feature (i) no longer causes problems, since we may apply 
the algorithm with a decreasing sequence of grid sizes, using the approximate solu- 
tion just found as a starting point for the next application. Algorithms of this 
form we shall call second-generation algorithms; they include Merrill's restart 
method [48,49] and Kuhn and MacKinnon' s sandwich method [41], discovered later but 
independently. 

The second method, introduced by Eaves [14], eliminates features (i) and (ii) 
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above. The grid size is automatically refined during the course of the algorithm. 
Such algorithms we will call third-generation, although they appeared at about the 
same time as Merrill’s. These methods are most easily explained after one is famil- 
iar with Merrill's algorithm; they are the subject of the next chapter. 

To motivate Merrill's algorithm, consider the following situation of Eaves’ 
first algorithm: 

Artificial labels based on c. 

In order to solve the right problem, we had to label vertices in C naturally. 
To obtain a starting point, artificial labels were necessary and had to be outside 
C. Hence the starting point was outside C. 

The artificial column d was needed for two reasons: 

(i) to obtain the required degree of freedom in the linear system; and 

(ii) to allow us to find a "deformed fixed point" (c+d) outside C and thus 
obtain a starting point. 

If we want to start within C, we must have artificial labels within C, in con- 
flict with our desire to have only natural labels in C. To solve the problem, we 
use two copies of C. 





C 



artificial 

labels 



C 



It is then natural to work in the (n+1 )-dimensional set [0,1] X C. Both reasons for 
d disappear. Since the dimension is increased by one, so is the number of vertices 
in a full-dimensional simplex, and we have the required degree of freedom. Also, we 
can start at a true fixed point of the artificial map of the bottom layer — it need 
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not be deformed by translation by d. 

It is clear why such algorithms are called sandwich methods — the two layers of 
the sandwich are {0} x C (artificial) and {l} x C (natural). (Eaves' first algo- 
rithm can be thought of as an "inefficient sandwich" where a hole is cut in a slice 
of bread for insertion of a piece of ham.) Using the same analogy, we can liken the 
third-generation algorithms to "club sandwiches": 



CO 



natural labels, 
finer grid 



A 



natural labels , 
coarse grid 






artificial 

labels 



VIII. 2 Preliminaries . We wish to find an approximate fixed point of a u.s.c. 

T1 

mapping F^: R -*■ R . Assume that satisfies the hypotheses of V.3.4: 

2.1 Assumption . There exist x^ £ R n , 6 > 0 , and y > 0 such that whenever 

x £ B(x°,y), f £ F 1 (x), and z £ B(x,6), (f-x) T (x°-z) > 0. 

Pick c £ R n arbitrarily, and let A be a non-singular n*n real matrix. 
Define F Q : R n -* R n ‘* by F Q (x) = {x + A(c-x)}. For any x = (x Q ,x i , . . . ,* n ) T € 

[0,1] xR n , let p(x) = (x^ , . . . ,x^)^ be its projection in R n . Then we can define 

F: [0,1] x R n + R n * by F(x) = x Q F 1 (p(x)) + (l-x 0 )F Q (p(x) ) . We call x a fixed 

T T 

point of F if F(x)=p(x). Note that (0,c ) is the only fixed point of F in 
{0}xR n and that any fixed point of F in {1} x R n projects into a fixed point of 



Let G be a triangulation of [0,1] x R n with C {0,1} xR n . Such a triangu- 
lation is called special. If a £ G, diam^a = sup{||p(x) - p(z)||^|x,z £ a} is the 



diameter of the projection of a, and mesh^G 



sup „ 0 diam a for q = 2, a 
*a£G q ^ 



The triangulation of R n whose simplices are (p(y ),.••, p(y D )) for 
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o = (y°,...,y n ) £ G n with o C (i)x R n is denoted G_^ , i = 0,1. 

a. 'b 

2.2 Examples . The restrictions and of and (as triangula- 

tions of R n+1 ) to [0,1] x R n are special triangulations. Let E denote the 
(n+l)x(n+l) matrix [v° ,ev\ . . . ,ev n ] ; then ) = {(Ey \...,Ey n )Ky \...,y n ) € 

'Xj % a. 

K^} and J^(e), defined similarly, are special triangulations. (K^e)^ = 

(K.(e)). = eK. and similarly for J, . We have mesh’ (K.. (e ) ) = mesh’ ( J. (e ) ) = e 
111 J 1 00 1 00 1 

and mesh^K^e ) ) = mesh^J^e)) = e/n. 

We will assume that c lies in an n-simplex of G^ or, equivalently, that 
T T 

(0,c ) lies in an n-simplex of G. (If not, c can be perturbed.) 

2.3 Definition . Let £ be a vector labelling for F with respect to G. If 
0 - (y \...,yk) £ G^ + ^, define the label matrix of o to be 



L 

a 



1 ... 1 

£(y' 1 ).--l(y k ) 



We say a is complete if there is a solution to 

L w = v^ , w > 0 , 
a — 



and very complete (v.c.) if there is a solution to 



(*) 



L a W = I, W > 0. (**) 

2 . 4 Lemma . Let F, G, and £ be as above and let a £ G n U G n+ \ Then there 
is a zero of £ in a iff a is complete. The only zero of £ in {0} x R n is 
(0,c T ) T , and the projection of any zero of £ in {l} x R n is a fixed point of a 
p.l. approximation to F^ with respect to G^. 

Proof . Similar to VII. 1.1. Q 



Note that £(y) depends only on c and A if y 6 {0} x R n and only on F^ 
if y £ { 1 } x R n . For this reason many authors define two labelling rules. Our 
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approach defines i throughout [0,1] xR n ; thus Merrill’s algorithm can be seen to 
trace zeroes of l from the artificial level to the natural level. 

VIII. 3 The Graph T . 

3.1 Definition . The nodes of T are v.c. (n+l)-simplices of G and n- 
simplices of G n lying in {0}xR n or (l}xR n . Two nodes are adjacent if one is 
a face of the other or if they share a v.c. facet. 

The following result shows that T has the right properties to be a basis for 
an algorithm. 

3 . 2 Lemma . 

(a) There is exactly one node of T lying in {0}xR n , namely, the n-simplex 

n T T 

Tq of G containing (0,c ) . 

(b) Any node of T lying in {1} x R n gives an approximate fixed point of F^. 

(c) If mesh^G <_ 6 and A is positive definite, and assumption 2.1 holds, 

T has only finitely many nodes. 

Proof . Parts (a) and (b) follow from 2.4; the fact that the simplex of G n 
T T 

containing (0,c ) is v.c. follows as in VI I. 3. 2. 

(c) We show that outside a compact region there are no v.c. simplices. We 
have x^, y, and 6 so that whenever x £ B(x^,y), f £ F^(x), and z £ B(x,<5), 
(f-x) T (x°-z) > 0. 

Let ||A|| 2 = max{ j|Ax[| 2 | ||x|| 2 = 1} and n = min{x T Ax j ||x || 2 = 1). Since 
{x £ R n |||xj| 2 = 1} is compact, the optima are attained and || A || 2 < °°, n > 0. 

Let y' = max{y + 5, j| A || 2 ( 6 + ||c-x° || 2 )/n } . Let C = [0 ,1] x B(x° ,y ’ ) C R n+1 . 

Then C is compact — we show that all complete simplices lie in C. Let a £ G have 

some vertex, y*, say, outside C; let z* = p(y*) and s = x^ - z*. We show 
T 

that s £(y) > 0 for all vertices y of a; hence a cannot be complete. 

Let y be a vertex of a with y Q = 1. Let z = p(y); we have £(y) = f - z 

for some f £ F^(z). Since z* £ B(x^,y+5) and z £ B(z* ,6), z £ B(x^,y). Hence 

T TO 

s £(y) = (f-z) (x -z*) > 0 by assumption 2.1. 

Let y be a vertex of a with y^ = 0 . Let z = p(y) — then £(y) = A(c-z) 
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and s T Jl(y) = (x°-z*) T A(c-z) = (x°-z*) T A(x°-z*) + (x°-z*) T A(c-x° ) + (x°-z*) T A( z*-z) . 
The first term above is positive; we claim that it dominates the other two terms. 
Indeed, if ||x^-z*|| 2 = v > y' , 

(x°-z*) T A(x°-z*) v 2 n > nvy’ >_ n v ||A || 2 ( 6 + ||c-x° || 2 ) /n . 

v ||A j| 2 5 is at least v||A(z*-z) || 2 and v||a|| 2 ||c-x°|| 2 is at least v||A(c-x°) || 2 . This 
T 

proves our claim; hence s i( y) >0 in this case also. 

Thus all complete simplices lie in C and the result is proved. Q 

Exactly as in VII. 3 we obtain that nodes of T have degree 1 or 2 according to 
whether they are n-simplices or (n+1 ) -simplices . 

3.3 Theorem . Each connected component of the finite graph T is either a 
simple circuit or a simple path each of whose endpoints is x^ or a v.c. n-simplex 
in {l}xR n . Q 

VI I I. 4 The Algorithm . 

4.1 Merrill's Basic Algorithm . Suppose we have F , <5, G, and l as in 
Section 2 with mesh^G <_ 6 . Let x Q be the simplex of G n containing (0,c T )^ 
and a the simplex of G n+1 with x Q as a facet. 

Step 0 . Form the lexicographic linear system showing x Q v.c. Let y + be 

the vertex of not in x^. Set m 1. 

Step 1. Calculate &(y + ) and introduce the column ( ^ ) into the current 

«.(y + ) 

basis , displacing a unique column ( ^ ) with y a vertex of a . 

«y"> .... n 

Step 2. If the facet x of opposite y lies in {1} x R , ST0P--T 

contains an approximate fixed point of F^. Otherwise, let a m+ j_ be the unique 

simplex of G sharing the facet x with a . Let y + be the new vertex of 

a , , set m m+1, and return to Step 1. 

m+1 

4.2 The Restart Algorithm . 

Step 0. Pick a sequence <5^ 4 0 with <_ 5 . For each k pick a special 
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triangulation with mesh’G^ = 6 V . Choose x° £ R n and set p «- 0 . 



(k) 



Step 1 . Apply algorithm 4.1 with ^(p) replacing G, replacing c 

some positive definite matrix to obtain an approximate fixed point x^ + \ 

p «- p+1 and return to Step 1. 

Of course, x^ must be perturbed if it does not lie in an n-simplex of 

(G ( P )V 

'Vi 

One typically takes G^^ = K^(a 6^) where 0*l<_a<_0-6. We choose a 
for "smooth" functions and large for highly nonlinear functions or mappings. 



and 

Set 



small 



VI 1 1. 5 An Example . Let n = 2 and F^x) 



10 0\,.2 l.T 

o o] (( 3’3 ) - x) ) ■ Clearl y> 



2 1 T 

has just one fixed point at (— ,— ) . Let us choose 



guess) and A = 



1 0 
0 1 



,2 

(y,y) (a perfect 



We will apply algorithm 4.1 using K^. 



We find x 
1 1 



, -1 0 1. / °\ M M\ 

= <y ,y >y > = M 0 )’! 1 J’l 1 ) )' 



The corresponding label matrix is 

-1 



this is the initial basis, with inverse 






0 

1 1 
0 -1 



Since the rows are lexicographically positive, x is v.c. The iterations now 
proceed as follows . 
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The course of the algorithm is shown in the following diagram. The dotted line 
is the locus of zeroes of l. 




This problem required 4 evaluations of a point of F^(x) and one evaluation of 

A(c-x), for a total of 5 linear programming pivot steps. In contrast, if we had 

taken A = l ) (so that the artificial map F_ coincides with F ), only 3 

\0 1 / 0 1 

function evaluations and 3 linear programming pivots would have been required. 

„ t / -1 0 1 2 v t,-1012 x 

Indeed, the simplices encountered are now = \y ,y ,y ,y = \y ,y >w ,y 

'-1012 10 
and a = (y ,w ,w ,y ); to see this, note that now &(^) = £(^) for any x. 

2 

Clearly, to obtain a v.c. simplex in {1} x R , at least three pivots are 
necessary (in general, we need n+1). The example illustrates this important point: 
to obtain fast convergence, it is not enough to make a good guess for c — we must 
also make a good guess for A. A should be close to the Jacobian of g if g is 
differentiable and F^x) = (x + g(x)}. In the early stages of applying algorithm 
4.2 little information about this Jacobian is available, and A must be positive 
definite to assure convergence. (Usually, A is taken to be the identity.) Later, 
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when fairly good approximate points have been generated, there is less concern that 
the algorithm will diverge — then we may take A to be an approximation to the 
Jacobian, even if it is not positive definite. The way to obtain such an approxi- 
mation without much additional work is the subject of exercise 7.5. The use of the 
Jacobian was first proposed in [55]; other references on its use are [20], [27], 

[66], [67] and [74]. The method for obtaining an approximate Jacobian appeared in 
[66] and [53]. 

VIII. 6 Integer Labelling . If = {f^} where f 1 is a continuous function, we 
can use Merrill's strategy together with integer labelling. We need a flexible 
artificial labelling to assure one and only one completely labelled simplex in 
{0}xR n . In fact, the following scheme achieves our goals. Define l as before. 
Then set the integer label of y to be 0 if l(y) < 0, and otherwise to be 
min{ j € N|A (y) = max^fc^y)}. 

To assure that the algorithm does not diverge, we must assume that every point 
outside a bounded region has a neighborhood of radius <5 in which some integer label 
is excluded. Since we want this condition to hold for any starting point c, we 
require the following. 

6.1 Assumption . There exist x° £ R n and 6 > 0 such that corresponding to 
every v > 0 there is some y > 0 with every x £ B(x^,y) satisfying: 
either (a) x - vu and f^z) £ z for all z £ B(x,6); or 

(b) for some i € N, x^ > x? + v and f}(z) - z^ < max^Cf^z) - z^) for 
all z £ B(x,6). 

Under this assumption algorithm 4.1, using integer labelling and a special 
triangulation G with mesh'G <_ 6, cannot diverge. For if y corresponds to 
v = 5 + ||x°-c|L o , any simplex outside [0 ,1] x B(x° ,y ) is missing label 0 (in case 
(a)) or label i (in case (b)). We omit the details. See [20] and [27] for related 
labelling rules. 
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VIII. 7 Exercises. 



7.1. Let L = 
a 



1 ... 1 

0 



c-y 



*c-y 



where a - (y^,...,y n ) is an n-simplex. Find L 



-1 



explicitly for a = k^(y^,ir) and o = j^(y^ ,tt ,s) . 

7.2 . Apply algorithm 4.1 exactly as in Section 5 but with F^(x) = 

r / 4 A 

{x + ( ) (c-x)}. Note: degeneracy occurs and lexicographic rules must be used. 

7.3 . Apply the integer labelling version of algorithm 4.1 to the example of 

Section 5 and with F, (x) = {x + ( ) (c-x)}. Note that in the latter case the 

1 \0 10/ 

final simplex does not contain the fixed point of F^. 

7.4 . State and prove a result giving the degree of approximation obtained when one 
uses the integer labelling of Section 6. 

7.5 . Imagine you are trying to find a fixed point of f, where f(x) = 

2 1 - 2 \ / 2 

2 3 -3 x + 3 J 9 which is equivalent to finding a zero of g with g(x) = 

V 3 - 5 / \3 



f(x)-x = Ax + b, Pi- 



ll 1 - 2 ' 

2 2-3 

\ 4 3 " 6 , 



Somehow you have found the 



simplex in 2K^ given by a = ((2,-4,0) T ,(2,-2, 0) T , (2,-2, 2) T , (4,-2, 2) T ) 
T 

k^( (2 ,-4 ,0) , (2,3,1)). The label matrix associated with a is 



” 1 1 1 1~ 




0 -X f 


0 2-20 


-1 


0 1 k 


-13-31 


with B = 


0 - - - — 
2 2 2 


-! 5 -7 1^ 




1 -i o I 

2 2 2 



Since the first column of B 






is nonnegative, we have a linear 



10 13 t 

approximate fixed point of f, given by — y + — y =(3,-3,l). Since f is 
linear, this is in fact a fixed point of f. However, much more information can be 
extracted from B 1 (which is available to us if a fixed-point algorithm gave a ) . 
First discard the first column of B -1 , then premultiply the rest by 
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0 111 



0 0 0 1 



° 1 'I 



R = 0 0 1 1 to get C = -1 — 0 . Now put the jth row of C into the 



-- 0 i 

2 2 



tt ( j ) T th position to get D 



Finally, multiply all entries in D 



"-3 0 1 

by the grid size 6=2 to get E = 0 2 -1 Now note that E = A 

_-2 1 0 _ 

Find out what operations on B give in turn C D \ and finally E ^ . Now 



apply these operations when B has the general form: 



g(y ) • • *g(y ) 



<y° , . . . ,y n ) = k^Cy 0 ,^), a simplex of 6K^. Conclude that the inverse of the label 
matrix can be used to find the inverse of an approximation to the Jacobian of g. 




CHAPTER IX: HOMOTOPY ALGORITHMS 



The third-generation algorithms we consider here are based on the work of Eaves 
[14] and Eaves and Saigal [16]. To motivate the methods, consider a continuous func- 

11 i 2 

tion f: S S . Instead of triangulating S , we triangulate S and omit one 

2 1 
simplex — S 2 is identified with S : 

v 




Label each vertex 0 or 1 according to whether its projection from v onto 
2 1 

S 2 = S has label 0 or 1. Then trace the natural path from the top 1-simp lex down 

to S 1 . The horizontal levels can be considered as successively finer triangulations 

of S 1 . These triangulations have meshes 1, -p, ^ above. 

A o 4 b b 

As another example, an alternative way to view the well-known bisection algo- 




Here the successive triangulations of S ^ have meshes 1 » ^ * * • • 

In both cases above, we used integer labelling for simplicity. However, if 
vector labelling is used, then the path links fixed points of p.l. approximations to 
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f with successively finer grids. The triangulation is used to deform f from a 
linear function to a p.l. approximation to f. 

IX. 1 Homotopies and Triangulations with Continuous Refinement of Grid Size . 

1.1 Definition. Let f 1 : R n -► R n be continuous functions for i = 0,1. Then 
h: [0,1] x R n -*■ R n is a homotopy from f^ to f^ if h is continuous and 
h(i,x T ) T = f 1 ( x ) for all x £ R n , i = 0,1. 

Note that the vector label i of Chapter VIII was a homotopy from A(c--) to 
a piecewise linear approximation to minus the identity. We traced zeroes of 

this function l. We would like our algorithms to trace fixed points of a homotopy 
from a given linear function to the function of interest. For computation these 
homotopies should be piecewise linear. 

Let f: R n -> R n and F = {f } : R n -> R n ". The constructions we will make are 

n n« 

motivated by F of this form but are also useful if F : R + R is merely u.s.c. 
We will note in brackets the differences for the latter case. 

The basic idea is to consider different levels on which are defined better and 
better approximations to f: 



level 


0 


f°, 


, a linear function 




level 


1 


f\ 


, a p.l. approx. 


to 


F with respect to a coarse 








triangulation 






level 


2 


f 2 ; 


> a p.l. approx. 


to 


F with respect to a finer 








triangulation 






level 


CO 


f 


itself 






:: if 


F * {f }, 


level 00 


is omitted.] 







It is convenient to put the levels in a geometric setting by considering level k 
as R n (k) = (2 k}*R n ; denote by R n (°°) the set {0} x R n . 

We now define the triangulations we use to construct the homotopies. 

1.2 Definition . Let 6^ V 0. Let G be a triangulation of (0,1] xR n 



satisfying 
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(i) if y £ G°, y Q = 2 _k for k= 0,1,2,...; 

(ii) for k = 0,l,2,... = {a £G n |aCR n (k)} triangulates R n (k); 

(iii) if ac(0,2^]xR n , diam’a <_ 5^_ for k = 0 ,1,2, . . . . 

Then G is called a triangulation with continuous refinement of grid size 
(abbreviated to "a triangulation with crogs"). 



1.3 Examples (n = 1). 




1 . 4 Lemma . If G is a triangulation with crogs and a £ G, then a c 
[2 k ^ # 2 k] x R n for some k >_ 0. 

Proof . Assume the contrary, i.e., that a has vertices y £ R D (p) and 
y f £ R n (^)» with p < k < l three integers. Then some convex combination z of 
y and y' lies in R n (k). Thus z £ t for some face t of a. This contradicts 
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the fact that z £ p for p some face of a' £ G^, from 1.2(ii). Q 

Using G we can construct a homotopy h from a linear function f^ to f. 

Let a: R + be a nondecreasing function, and let f^(x) = A(c-x) + x. 

1.5 Definition . Define F T : (0,1] xR n -* R n “ by F’ (x) = {f^(p(x))} if 

-log 2 x Q <_ a(||p(x) - c || 2 ) and F f (x) = F(p(x)) otherwise. Restricted to (0,l]xR n , 

T T 

h is a piecewise linear approximation to F' with respect to G. Set h(0,x ) = 

f (x) ; thus h: [0,l]xR n -► R n . [If F i { f } , omit the last sentence; we have 
h: (0,1] xR% R n .] 

We usually set a identically to zero; then we label according to f^ on 
R n (0) and F on all other levels. 

1.6 Lemma , h is a homotopy from f to f°. [1.6 is meaningless if F / 

{f }- ] 

T T T T 0 . o 

Proof . Clearly, h(0,x ) = f(x) and h(l,x ) = f (x) since f is linear 

and agrees with its p.l. approximation. Also, h is continuous on (0,1] x R n 

because it is a p.l. approximation. So let x £ R n (°°) and e > 0 be given. The 

continuity of f gives 6 > 0 such that ||f(w) - f(p(x))|| 2 < e whenever 

w £ B(p(x),6). There exists k > a(||p(x) - c|| + 6) so that 6^ <_ 6/2. Let M = 

[0,2 ^ ■*"] x B(p(x) ,6/2) and z £ M. We claim that ||h(z) - h(x) || <_z. 

If z Q = 0, h(z) = f(p(z)) and the claim is true. If z Q > 0, z lies in 

a for a = <y _1 ,...,y n ) £ G. By 1.4, y^ <_ 2 k for all i and hence 1.2(iii) 

gives diam'a £ 6^ <_ 6/2. Thus UpCy 1 ) - p(x) || 2 <_ 5/2 + 6/2 = 6 and 

||h (y 1 ) - h(x)|| 2 = ||f (pCy 1 ) ) - f(p(x))|| 2 for all i. The claim follows. [~~| 

1.7 Definition . We call x £ [0,1] xR n a fixed point of h if h(x) = p(x). 

Thus x £ R n (0) is a fixed point of h iff p(x) = c or x = (l,c T ) T when A 

is nonsingular. Also, x £ R n (°°) is a fixed point of h iff p(x) is a fixed 
point of F. [The last statement is meaningless if F 4 { f } . ] 

A finite algorithm cannot obtain a fixed point of h in R n (°°). The following 



result is therefore of interest. 
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1 . 8 Lemma . Let {x } be a sequence of fixed points of h tending to 
x £ R n (°°). Then p(x) is a fixed point of F. 

Proof . If F = {f}, we have h: [0,1] xR n R n continuous, so h(x) = 
k k 

lim h(x ) = lim p(x ) = p(x). [Otherwise, the conclusion follows from 1.2(iii) and 
the method of proof of Kakutani's theorem.] | [ 



IX. 2 Complete Simplices . 

2.1 Definition . Define l: [0,1] x R n R n by £(x) = h(x) - p(x). [If 
F i {f}, l : (0,1] x R n + R n is defined similarly.] If a £ G n U G n+1 , the label 
matrix of a is 



L 

a 



i( y _1 )...£(y k ) 



where a = (y ^,...,y^). 



We call a complete (or very complete) iff there is a feasible solution to L^w = 

v°, w 0 (L W = I, W > 0). 

2.2 Lemma . Let cr £ G n U G n+ '*". Then there is a fixed point of h in a iff 

. T T — 

a is complete. If a £ G Q is complete, (l,c ) £ a. 

Proof. Trivial. 



IX. 3 The Graph T . 

3.1 Definition . The nodes of T are v.c. (n+l)-simplices and v.c. n- 
simplices of G Q , with two such adjacent if one is a face of the other or if they 
share a v.c. facet. 

T T 

Assume that (l,c ) lies in an n-simplex of G Q . Then it follows that 

Tq is the unique node of Y that is an n-simplex, and T contains an infinite 
simple path with as its endpoint. The algorithm will follow this path. We want 

to assure that the path eventually penetrates every level. 

Assume that F satisfies assumption 2.1 of Chapter VIII with x° , u and 
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(5 > 0. Suppose A is positive definite. Then as in VIII. 3.2 no complete simplex 
of diameter at most 6 can lie outside f = (0,l] x B(x ,y T ) for some u ' . 

Thus if G is a triangulation with crogs with meshes anc * ^0 — ^ ’ a11 

nodes of F lie in C. By appropriate use of the function a in 1.5 we obtain a 
stronger result. 

3.2 Lemma. Under the conditions above, with 5 q <_ 6, all nodes of T lie in 

C^,. Furthermore, if a is an increasing function, then even with 6^> 6 , all 

nodes of T lie in C for some y" > 0. 

V" 

Proof. If 6^ < 6, the proof follows that of VIII. 3.2. Otherwise, pick k 

so that 5^ <_ 5 and y" > u' so that a(y" - ||x°-c|| 2 - <5 Q ) >_ k+1. Let a be a 

*~]c n 

simplex with some vertex outside C^„. If o C (0,2 ] X R , o is not complete 

as in VIII. 3. 2. If a C [2~ k ,l] x R n , x £ a => -log^ < k+1 and hence each 

vertex y of a is labelled A(c-p(y)). Thus if a is complete, it contains 
(X,c T ) T for some A, which contradicts a £ C . Q 

Note that the same idea can be used in Merrill's algorithm to assure conver- 
gence even if 6^ > 6. 



IX. 4 The Algorithm . We trace the infinite path of T from t q until a suffi- 
ciently accurate approximate fixed point is found. By 3.2 any infinite sequence of 
distinct, very-complete simplices must leave [2 ,1] x R for any k; by 2.2 we 
obtain a fixed point of h in R n (k) for every k; and, by 1.8, any cluster point 
of these fixed points is a fixed point of F. 

Step 0. Let be the unique (n+l)-simplex of G with as a facet. Let 

y + be the vertex of a outside R n (0). Form the lexicographic linear system 
showing Tq very complete. Set m +- 1. 

Step 1. Calculate &(y + ). Introduce the column ( ) into the current 

«.(y ) -l 

basis of the lexicographic linear system. Assume that the column _^) drops 

from the basis , with y a vertex of a . 

Step 2. Let be the very-complete facet of opposite y . If desired, 

calculate the fixed point of h in i m and terminate if it is sufficiently 
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accurate. Let a be the simplex of G sharing the facet x with c . Let 
m+1 m m 

y + be the new vertex, set m «- m+1, and return to Step 1. 

Generally, the fixed point of h in x^ is calculated only when x^ C R n (k) 
for some k. 

IX. 5 A Comparison with Merrill’s Algorithm . The advantages of a homotopy algorithm 
over Merrill's algorithm are: 

1) The function is continuously deformed rather than set to a constant function for 
each restart; and 

2) The simplices generated can move from fine triangulations back to coarse ones, 
thus permitting a move to a more promising region in large simplices. 




To eliminate advantage (1), we can use A(c-x) as the artificial label, with 
A an approximation to the Jacobian of f minus the identity (see VIII. 5 and 
VIII. 7. 5). The second supposed "advantage" has been discouraging in computational 
experience — generally, Merrill's algorithm has the advantage when the homotopy 
algorithms "yoyo”. 

The advantage of Merrill's algorithm is its greater flexibility. The grid size 
can be reduced by any factor, as compared to i for presently known triangulations 
with crogs. Also, one can make more use of the approximate inverse Jacobian in 
taking quasi-Newton steps with the option to restart the algorithm if convergence 
is poor. Of course, one can use a mixed strategy, and restart a homotopy algorithm 
as often as desired. Saigal, Solow, and Wolsey have a flexible algorithm [55] 
incorporating advantage (1) of homotopy algorithms with the flexibility of Merrill's 
approach . 




CHAPTER X: TRIANGULATIONS WITH CONTINUOUS REFINEMENT OF GRID SIZE 



X.l Discussion . Triangulations with crogs were first introduced by Eaves [14] and 
Eaves and Saigal [16]. These triangulations, known as K^, and Kg, differ from 
our triangulations of the same name. Diagrams of Eaves’ and are in IX. 1.3. 

Recall the following: 

1.1 Definition . G is a triangulation with crogs if for some sequence 6^ * 0, 

(i) G triangulates (0,1] xR n and G° C R n (k); 

(ii) = {a £ G n |a c R n (k)} triangulates R n (k) for k = 0,1,..,; and 

■“k n t 

(iii) o C (0 ,2 ] x R implies mesh 2 a <_ 6^ for k = 0 ,1 , . , . . 

A natural way to construct such triangulations is to choose G^ to be 
or with and then join up the G^’s suitably. The first ques- 

tion is to choose the e^. Nobody knows yet how to construct triangulations with 
crogs where e < \ (except for, say, n = 1, when any ratio is possible.) 

However, even if one could construct such triangulations, one would incur some verti- 
cal paths between levels of lengths greater than n+1. The importance of this Result 
is that with an affine function (or with a continuously differentiable function and 
sufficiently small mesh) a homotopy algorithm follows a vertical path. 

1 . 2 Lemma . Let G be a triangulation with crogs so that mesh 2 G^_ + ^/mesh 2 G^ < 
~ for some k. Then there exists x £ R n such that [2 ^ \ 2 ^] x {x} intersects 
more than n+1 simplices of G. 

Proof . Let a be a simplex of G^ with diam^a > 2 mesh^G^^^, and let x 
and z be two points of a with |(x-z ||^ > 2 mesh^G^^^. Let y £ R n (k+1) be the 
vertex of a simplex p of G with a as a facet. Then either ||p ( x ) - p(y) |l 2 or 

||p(z) - p(y) || 2 is greater than mesh^G^^^; assume the first. Now the sequence of 
— k “■ 1 -k 

simplices met by [2 ,2 ] x {x} starts with a simplex with n+1 vertices in 

R n (k+1) with none equal to y and ends with a simplex with y the only vertex 
in R n (k+1). Since each pivot transfers at most one vertex from R n (k+1) to R n (k), 
at least n+2 simplices are met. Q 



It would be useful to find a strengthening of 1.2 that related the maximum 




